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Abstract 

The condensation of elementary quanta and their macroscopic occupation of the same quan- 
tum state, say k = in some reference frame S, is the essential ingredient of the degenerate 
vacuum of present-day elementary particle physics. This represents a sort of 'quantum ether' 
which characterizes the physically realized form of relativity and could play the role of pre- 
ferred reference frame in a modern re-formulation of the Lorentzian approach. In spite of this, 
the so called 'null results' of the classical ether-drift experiments, traditionally interpreted as 
confirmations of Special Relativity, have so deeply influenced scientific thought as to prevent 
a real critical discussion on the reasons underlying its alleged supremacy. In this paper, we 
argue that this traditional null interpretation is far from obvious. In fact, by using Lorentz 
transformations to connect the Earth's frame to S, the small observed effects point to an 
average Earth's velocity of about 300 km/s, as in most cosmic motions. A common feature 
is the irregular behaviour of the data. While this has motivated, so far, their standard inter- 
pretation as instrumental artifacts, our new re-analysis of the very accurate Joos experiment 
gives clear indications for the type of Earth's motion associated with the CMB anisotropy and 
leaves little space for this traditional interpretation. The new explanation requires instead a 
view of the vacuum as a stochastic medium, similar to a fluid in a turbulent state of motion, in 
agreement with basic foundational aspects of both quantum physics and relativity. The over- 
all consistency of this picture with the present experiments with vacuum optical resonators 
and the need for a new generation of dedicated ether-drift experiments are also emphasized. 
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1. Introduction 



An analysis of the ether-drift experiments, starting from the original Michelson-Morley exper- 
iment of 1887, should be suitably framed within a general discussion of the basic differences 
between Einstein's Special Relativity [1] and the Lorentzian point of view [11 [3l |3]. There 
is no doubt that the former interpretation is today widely accepted. However, in spite of 
the deep conceptual differences, it is not obvious how to distinguish experimentally between 
the two formulations. This type of conclusion was, for instance, already clearly expressed by 
Ehrenfest in his lecture 'On the crisis of the light ether hypothesis' (Leyden, December 1912) 
as follows: "So, we see that the ether-less theory of Einstein demands exactly the same here 
as the ether theory of Lorentz. It is, in fact, because of this circumstance, that according to 
Einstein's theory an observer must observe exactly the same contractions, changes of rate, 
etc. in the measuring rods, clocks, etc. moving with respect to him as in the Lorentzian 
theory. And let it be said here right away and in all generality. As a matter of principle, 
there is no experimentum crucis between the two theories". This can be understood since, 
independently of all interpretative aspects, the basic quantitative ingredients, namely Lorentz 
transformations, are the same in both formulations. Their validity will be assumed in the 
following to discuss the possible existence of a preferred reference frame. 

For a modern presentation of the Lorentzian philosophy one can then refer to Bell [5l 
[6l [?]• In this alternative approach, differently from the usual derivations, one starts from 
physical modifications of matter (namely Larmor's time dilation and Lorentz-Fitzgerald length 
contraction in the direction of motion) to deduce Lorentz transformations. In this way, due to 
the fundamental group properties, the relation between two observers S' and S" , individually 
related to the preferred frame S by Lorentz transformations with dimensionless parameters 
/3' = v' /c and /3" = v" /c, is also a Lorentz transformation with relative velocity parameter 
/3rei fixed by the relativistic composition rule 



(for simplicity we restrict to the case of one-dimensional motion). This produces a substan- 
tial quantitative equivalence with Einstein's formulation for most standard experimental tests 
where one just compares the relative measurements of a pair of observers. Hence the impor- 
tance of the ether-drift experiments where one attempts to measure an absolute velocity. 

At the same time, if the velocity of light propagating in the various interferometers 
coincides with the basic parameter c entering Lorentz transformations, relativistic effects 
conspire to make undetectable the individual /?', /?",... This means that a null result of the 
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ether-drift experiments should not be automatically interpreted as a confirmation of Special 
Relativity. As stressed by Ehrenfest, the motion with respect to S might remain unobservable, 
yet one could interpret relativity ' a la Lorentz'. This could be crucial, for instance, to reconcile 
faster-than-light signals with causality [8j and thus provide a different view of the apparent 
non-local aspects of the quantum theory j9j. 

However, to a closer look, is it really impossible to detect the motion with respect to S? 
This possibility, which was implicit in Lorentz' words [5] "...it seems natural not to assume 
at starting that it can never make any difference whether a body moves through the ether or 
not..", may induce one to re-analyze the classical ether-drift experiments. Let us first give 
some general theoretical arguments that could motivate this apparently startling idea. 

A possible observation is that Lorentz symmetry might not be an exact symmetry. In 
this case, one could conceivably detect the effects of absolute motion. For instance Lorentz 
symmetry could represent an 'emergent' phenomenon and thus reflect the existence of some 
underlying form of ether. This is an interesting conceptual possibility which, in many different 
forms, objectively reflects the fast growing interest of part of the physics community, a partial 
list including i) the idea of the vacuum as a quantum liquid |10[ [TT] (which can explain in a 
natural way the huge difference between the typical vacuum-energy scales of modern particle 
physics and the cosmological term needed in Einstein's equations to fit the observations) 
ii) the idea of Lorentz symmetry as associated with an infrared fixed point \12\ [T3] in non- 
symmetric quantum field theories iii) the quantum-gravity literature which, by starting from 
the original concept [14j of 'space-time foam', explicitly models the vacuum as a turbulent 
fluid \15\ [TBI [T7] iv) the idea of deformations of Lorentz symmetry in a theoretical scheme 
('Doubly Special Relativity') [181 13 ED] where besides an invariant speed there is also an 
invariant length associated with the Planck scale v) the representation of relativistic particle 
propagation from the superposition, at very short time scales, of non-relativistic particle paths 
with different Newtonian mass |21| . 

Here, however, we shall adopt a different perspective and concentrate our analysis on 
a peculiar aspect of today's quantum field theories: the representation of the vacuum as a 
'condensate' of elementary quanta. These condense because their trivially empty vacuum is 
a meta-stable state and not the true ground state of the theory. In the physically relevant 
case of the Standard Model of electroweak interactions, this situation can be summarized 
by saying [22] that "What we experience as empty space is nothing but the configuration of 
the Higgs field that has the lowest possible energy. If we move from field jargon to particle 
jargon, this means that empty space is actually filled with Higgs particles. They have Bose 
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condensed". The explicit translation from field jargon to particle jargon, with the substantial 
equivalence between the effective potential of quantum field theory and the energy density of 
a dilute particle condensate, can be found for instance in ref . j23] . 

The trivial empty vacuum will eventually be re-established by heating the system above 
a critical temperature T = Tc where the condensate 'evaporates'. This temperature in the 
Standard Model is so high that one can safely approximate the ordinary vacuum as a zero- 
temperature system (think of ^He at a temperature 10~^^ °K). This observation allows one 
to view the physical vacuum as a superfiuid medium pi)] where bodies can flow without any 
apparent friction, consistently with the experimental results. Clearly, this form of quantum 
vacuum is not the kind of ether imagined by Lorentz. However, if possible, this modern view 
of the vacuum state is even more different from the empty space-time of Special Relativity 
that Einstein had in mind in 1905. Therefore, one might ask [24j if Bose condensation, i.e. the 
macroscopic occupation of the same quantum state, say k = in some reference frame S, can 
represent the operative construction of a 'quantum ether'. This characterizes the physically 
realized form of relativity and could play the role of the preferred reference frame in a modern 
Lorentzian approach. 

Usually this possibility is not considered with the motivation, perhaps, that the average 
properties of the condensed phase are summarized into a single quantity which transforms 
as a world scalar under the Lorentz group, for instance, in the Standard Model, the vacuum 
expectation value (<&) of the Higgs field. However, this does not imply that the vacuum state 
itself has to be Lorentz invariant. Namely, Lorentz transformation operators C/', [/",.. might 
transform non trivially the reference vacuum state |^'('^)) (appropriate to an observer at rest 
in S) into j^*'), |^'")v (appropriate to moving observers 5', 5"',.-) and still, for any Lorentz- 
invariant operator G, one would find 

{G)^(o) = {G)-q,i = {G)-^" = .. (2) 

Here, we are assuming the existence of a suitable operatorial representation of the Poincare 
algebra for the quantum theory in terms of 10 generators Pa, M^^is ( a ,/3=0, 1, 2, 3) where 
Pa are the 4 generators of the space-time translations and Mai3 = —Mpa are the 6 generators 
of the Lorentzian rotations with commutation relations 

[Po.,Pp]=^ (3) 

[Map, P^] = Vl^jPa - Va^Pp (4) 
[Map, M^s] = rja-yMps + '^psMa-^ - Tlp-yMaS - VaSMp^ (5) 
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where r]ai3 = diag(l, —1, —1, —1). 

With these premises, the possibility of a Lorentz-non- invariant vacuum state was addressed 
in refs.|251[26] by considering two basically different approaches. In the first description, as in 
the axiomatic approach to quantum field theory [27], one could describe the physical vacuum 
as an eigenstate of the energy-momentum vector. This physical vacuum state 

1^(0)) 

would 

maintain both zero momentum and zero angular momentum, i.e. (i,j = l,2,3) 

A|*^°^> = Mij|*^°^) = (6) 

but, at the same time, be characterized by a non- vanishing energy 

Pol^(o)) = So|^^°^) (7) 

This vacuum energy might have different explanations. Here, we shall limit ourselves to explor- 
ing the physical implications of its existence by just observing that, in interacting quantum 
field theories, there is no known way to ensure consistently the condition Eq = without 
imposing an unbroken super symmetry, which is not phenomenologically acceptable. In this 
framework, by using the Poincare algebra of the boost and energy-momentum operators, one 
then deduces that the physical vacuum cannot be a Lorentz- invariant state and that, in any 
moving frame, there should be a non-zero vacuum spatial momentum (Pj)^' 7^ along the 
direction of motion. In this way, for a moving observer 5' the physical vacuum would look like 
some kind of ethereal medium for which, in general, one can introduce a momentum density 
{Woi)q,' through the relation (i=l,2,3) 

(Pi)*' = J d?x {Wor)^' ^ (8) 

On the other hand, there is an alternative approach where one tends to consider the vacuum 
energy Eq as a spurious concept and only concentrate on an energy-momentum tensor of the 
following form [28l [29] 

{W^l,)^(0) = Pv VfMU (9) 

{pv being a space-time independent constant). In this case, one is driven to completely 
different conclusions since, by introducing the Lorentz transformation matrices to any 
moving frame S", defining {Wf^u)^' through the relation 

{W^,h> = A%AP, (VF.p)*(o) (10) 

^ We ignore here the problem of vacuum degeneracy by assuming that any overlapping among equivalent 
vacua vanishes in the infinite- volume limit of quantum field theory (see e.g. S. Weinberg, The Quantum Theory 
of Fields, Cambridge University press, Vol.11, pp. 163-f67). 
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and using Eq.Q, it follows that the expectation value of Woi in any boosted vacuum state 
1^'') vanishes, just as it vanishes in |^'^'^^), i.e. 

j (fx {Woi)^, = {P^)^> = (11) 

As discussed in ref . [25] . both alternatives have their own good motivations and it is not so 
obvious how to decide between Eq.(l8|) and Eq. ijlip on purely theoretical grounds. For instance, 
in a second-quantized formalism, single-particle energies -E'i(p) are defined as the energies of 
the corresponding one-particle states |p) minus the energy of the zero-particle, vacuum state. 
If Eq is considered a spurious concept, Ei{p) will also become an ill-defined quantity. At a 
deeper level, one should also realize that in an approach based only on Eq.Q the properties 
of |^'('^)) under a Lorentz transformation are not well defined. In fact, a transformed vacuum 
state 1^'') is obtained, for instance, by acting on I'l'*^''^) with the boost generator Mqi- Once 
is considered an eigenstate of the energy-momentum operator, one can definitely show 
that, for / 0, 1^-') and |^'(°)) differ non-trivially. On the other hand, if £"0 = there 
are only two alternatives: either Afoil^'^"-') = 0, so that 1^') = l^-e^)), or Moi|*(°)) is a state 
vector proportional to l^''^'^-'), so that |\I'') and l^'^'^^) differ by a phase factor. 

Therefore, if the structure in Eq.Q were really equivalent to the exact Lorentz invariance 
of the vacuum, it should be possible to show similar results, for instance that such a |\I'(^^) 
state can remain invariant under a boost, i.e. be an eigenstate of 

Moi = -ij d^x {xMoo - xoWoi) (12) 

with zero eigenvalue. As far as we can see, there is no way to obtain such a result by 
just starting from Eq.Q (this only amounts to the weaker condition (Moi),j,(o) = 0). Thus, 
independently of the finiteness of Eq^ it should not come as a surprise that one can run into 
contradictory statements once |\I'('^)) is instead characterized by means of Eqs.® — d?]). For 
these reasons, it is not obvious that the local relations ([9]) represent a more fundamental 
approach to the vacuum. 

Alternatively, one could argue that a satisfactory solution of the vacuum-energy prob- 
lem lies definitely beyond flat space. A non-zero p^,, in fact, should induce a cosmological 
term in Einstein's field equations and a non-vanishing space-time curvature which anyhow 
dynamically breaks global Lorentz symmetry. Nevertheless, in our opinion, in the absence of 
a consistent quantum theory of gravity, physical models of the vacuum in flat space can be 
useful to clarify a crucial point that, so far, remains obscure: the huge renormalization effect 
which is seen when comparing the typical vacuum-energy scales of modern particle physics 
with the experimental value of the cosmological term needed in Einstein's equations to fit 
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the observations. For instance, as anticipated, the picture of the vacuum as a superfluid 
can explain in a natural way why there might be no non-trivial macroscopic curvature in 
the equilibrium state where any liquid is self-sustaining [lOj. In any liquid, in fact, curva- 
ture requires deviations from the equilibrium state. The same happens for a crystal at zero 
temperature where all lattice distortions vanish and electrons can propagate freely as in a 
perfect vacuum. In such representations of the lowest energy state, where large condensation 
energies (of the liquid and of the crystal) play no observable role, one can intuitively under- 
stand why curvature effects can be orders of magnitude smaller than those naively expected 
by solving Einstein's equations with the full {W^u)^{o) as a cosmological term. In this per- 
spective, 'emergent-gravity' approaches \30 \ \31 \ [32]. where gravity somehow arises from the 
same physical flat-space vacuum, may become natural and, to find the effective form for 
the cosmological term to be inserted in Einstein's field equations, we are lead to sharpen our 
understanding of the vacuum structure and of its excitation mechanisms by starting from the 
physical picture of a superfluid medium. To decide between Eqs.([8]) and (fTTj) . one could then 
work out the possible observable consequences and check experimentally the existence of a 
fundamental energy-momentum flow. 



2. Vacuum energy- momentum flow as an ether drift 

To explore the idea of a non-zero vacuum energy-momentum flow, one can adopt a phe- 
nomenological model [23] where the physical vacuum is described as a relativistic fluid |34] . 
In this representation, a non-zero (Woi)*' gives rise to a tiny heat flow and an effective ther- 
mal gradient in a moving frame S' . This would represent a fundamental perturbation which, 
if present, is likely too small to be detectable in most experimental conditions by standard 
calorimetric devices. However, it could eventually be detected through very accurate ether- 
drift experiments performed in forms of matter that react by producing convective currents 
in the presence of arbitrarily small thermal gradients, i.e. in gaseous systems. 

To better explain this possibility, let us first recall that in the modern version of these 
experiments one looks for a possible anisotropy of the two-way velocity of light through the 
relative frequency shift Az/(0) of two orthogonal optical cavities [35\ 136] . Their frequency 



^ In this sense, by exploring emergent-gravity approaches based on an underlying superfluid medium, one 
is taking seriously Feynman's indication : "...the first thing we should understand is how to formulate gravity 
so that it doesn't interact with the energy in the vacuum" ,33 • 
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is proportional to the two-way velocity of light Cy{9) within the cavity through an integer 
number m, which fixes the cavity mode, and the length of the cavity L(9) as measured in 
the laboratory. In principle, by filling the resonating cavities with some gaseous medium, the 
existence of a vacuum energy-momentum flow could produce two basically different effects: 

a) modifications of the solid parts of the apparatus. These can change the cavity length 
upon active rotations of the apparatus or under the Earth's rotation. 

b) convective currents of the gas molecules inside the optical cavities. These can produce 
an anisotropy of the two-way velocity of light. In this sense, the reference frame S' where the 
solid container of the gas is at rest would not define a true state of rest. 

Now, an anisotropy of the cavity length, in the laboratory frame, would amount to an 
anisotropy of the basic atomic parameters, a possibility which is severely limited experimen- 
tally. In fact, in the most recent versions of the original Hughes-Drever experiment |37[ I38|. 
where one measures the atomic energy levels as a function of their orientation with respect to 
the fixed stars, possible deviations from isotropy have been found below the IQ-^o level [39] . 
This is incomparably smaller than any other effect on the velocity of light that we are going 
to discuss. Therefore, mechanism a), if present, is completely negligible and, from now on, we 
shall assume L{9) = L =constant. In this way, one re-obtains the standard relation adopted 
in the analysis of the experiments 

AuP'-y^e) ^ c^{7r/2 + 9)-c^{d) _ Ace ^^^^ 

Uq C C 

where is the reference frequency of the two optical resonators and the suffix "phys" indicates 
a hypothetical physical part of the frequency shift after subtraction of all spurious effects. 

Let us now estimate the possible effects of mechanism b) by first recalling that rigorous 
treatments of light propagation in dielectric media are based on the extinction theory |4Uj . 
This was originally formulated for continuous media where the inter-particle distance is smaller 
than the light wavelength. In the opposite case of an isotropic, dilute random medium [51] 
as a gas, it is relatively easy to compute the scattered wave in the forward direction and 
obtain the refractive index. However, the presence of convective currents would produce an 
anisotropy of the velocity of refracted light. 

To derive the relevant relations, let us introduce from scratch the refractive index M of the 
gas. By assuming isotropy, the time t spent by refracted light to cover some given distance 
L within the medium is t = ML/c. This can be expressed as the sum of to = L/c and 
ti = (J\f — ^)L/c where to is the same time as in the vacuum and ti represents the additional, 
average time by which refracted light is slowed down by the presence of matter. If there are 
convective currents, due to the motion of the laboratory with respect to a preferred reference 
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frame E, then ti will be different in different directions, and there will be an anisotropy of 
the velocity of light proportional to {M — 1). In fact, let us consider light propagating in a 
2-dimensional plane and express ti as 



h = -f{M,9,l3) 
c 



(15) 



with f3 = V/c, V being (the projection on the considered plane of) the relevant velocity with 
respect to S where the isotropic form 



(16) 



is assumed. By expanding around M = 1 where, whatever /3, / vanishes by definition, one 
finds for gaseous systems (where A/" — 1 ^ 1) the universal trend 



/(AA,0,/3)~(AA-l)F(e,/3) 



with 



Fie,f3)^{df/dM)W=i 

and F{9,0) = 1. Therefore, by introducing the one-way velocity of light 



one gets 



c^{M, /3) ~ - [1 - (AA - 1) {F{e, (3) - 1)] 



(17) 
(18) 

(19) 
(20) 



Analogous relations hold for the two-way velocity Cy{M,0, f]) 



2 c^(AA,g,/3)c^(AA,7r + g,/3) 
c^{M,e,l3) + c^{M,7T + e,l3) ^ 'M 



21 



A more explicit expression can be obtained by exploring some general properties of the func- 
tion F{6,/3). By expanding in powers of /3 



Fie,p)-i = /3F,{e) + (3^F2{e) + ... 



(22) 



and taking into account that, by the very definition of two-way velocity, c^(A^, ^,/3) = 
c^{Af,0,-f3), it follows that Fi{e) = -Fi{tt + 6). Therefore, to 0{f3'^), we get the general 
structure [26] 



c^(AA,0,/3)~- 



n=0 



(23) 
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in which we have expressed the combination F2{9) + F2{Tr + 9) as an infinite expansion of 
even-order Legendre polynomials with unknown coefficients = which depend on the 
characteristics of the induced convective motion of the gas molecules inside the cavities. 

Eq. (j23p . in principle, is exact to the given accuracy but it is of limited utility if one wants 
to compare with real experiments. In fact, it would require the complete control of all possible 
mechanisms that can produce the gas convective currents by starting from scratch with the 
macroscopic Earth's motion in the physical vacuum. This general structure can, however, be 
compared with the particular form (see Eq. ljlOSp of the Appendix) obtained by using Lorentz 
transformations to connect S' to the preferred frame 

c^(AA,0,/3) --^[1-/32 (M-1){A + Bsm'^e)] (24) 

with A = 2 and B = —1 which corresponds to setting ^0 = 4/3, C2 = 2/3 and all C2n = 
for n > 1 in Eq. (j23|) . Eq.(|24|) represents a definite realization of the general structure in (|23p 
and a particular case of the Robertson-Mansouri-Sexl (RMS) scheme [4:2\ |33j for anisotropy 
parameter \B\ = M — 1 (see the Appendix). In this sense, it provides a partial answer to the 
problems posed by our limited knowledge of the electromagnetic properties of gaseous systems 
and will be adopted in the following as a tentative model for the two-way velocity of light 

Summarizing: in this scheme, the theoretical estimate for a possible anisotropy of the 
two-way velocity of light is 

(AAgas - 1) ^ (26) 

Then, by assuming the typical velocity of most Earth's cosmic motions V ~ 300 km/s, one 
would expect < 10~^ for experiments performed in air at atmospheric pressure, where 





Theor 


"Ace" 


Theor 


. ^0 . 


gas 


c 


gas 



■^One conceptual detail concerns the gas refractive index whose reported values are experimentally measured 
on the Earth by two-way measurements. For instance for air, the most precise determinations are at the level 
10"^, say A/kir ~ 1.0002926.. at STP (Standard Temperature and Pressure). By assuming a non-zero anisotropy 
in the Earth's frame, one should interpret the isotropic value c/A/lir as an angular average of Ea. (|24p . i.e. 

^ {c,(A4i„e,/3)>« = [1 - f (A/^ir - 1)/3'1 (25) 

From this relation, one can determine in principle the unknown value A/lir = A/'(E) (as if the gas were at 
rest in E), in terms of the experimentally known quantity A/kir = M {Earth) and of V . In practice, for the 
standard velocity values involved in most cosmic motions, say V ~ 300 km/s, the difference between A/'(E) 
and Af (Earth) is at the level 10~® and thus completely negligible. The same holds true for the other gaseous 
systems at STP (say nitrogen, carbon dioxide, helium,..) for which the present experimental accuracy in the 
refractive index is, at best, at the level 10~^. Finally, the isotropic two-way speed of light is better determined 
in the low-pressure limit where {J\f — 1) — >■ 0. In the same limit, for any given value of V, the approximation 
A/'(E) — J\f (Earth) becomes better and better. 
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M ~ 1.00029, or < 10~^^ for experiments performed in helium at atmospheric pressure, 
where J\f ~ 1.000035. Therefore these potential effects are much larger than those possibly 
associated with vacuum cavities. In fact, from experiments one finds [44j — [50] 



"Ai/' 


EXP 


"Acs' 


EXP 


. ^0 . 


vacuum 


c 


vacuum 



~ 10-15 (27) 



or smaller and thus completely negligible when compared with those of Eq. (j26p . 

On the other hand, if one were considering light propagation in a strongly bound system, 
such as a solid or liquid transparent medium, the small energy flow generated by the motion 
with respect to the vacuum condensate should mainly dissipate by heat conduction with no 
appreciable particle flow and no light anisotropy in the rest frame of the container of the 
medium. This conclusion is in agreement with the experiments [71 [51] that seem to indicate 
the existence of two regimes. A former region of gaseous systems where A/" ~ 1 and there are 
small residuals which are roughly consistent with Eq. (|26|) . A latter region where the difference 
of M from unity is substantial, (e.g. A/" ~ 1.5 as with perspex in the experiment by Shamir 
and Fox [52j), where light propagation is seen isotropic in the rest frame of the medium (i.e. 
in the Earth's frame). Although it would be difficult to describe in a fully quantitative way 
the transition between the two regimes, some simple arguments can be given along the lines 
suggested by de Abreu and Guerra (see pages 165-170 of ref. [53] ) . 

For this reason, it was proposed in refs.[71 [25l |26] that one should design a new class of 
dedicated experiments in gaseous systems. Such a type of 'non-vacuum' experiment would be 
along the lines of ref.|54j where just the use of optical cavities filled with different materials 
was considered as a useful complementary tool to study deviations from exact Lorentz invari- 
ance. In the meantime, due to the heuristic nature of our approach, and to further motivate 
this new series of experiments, one could try to obtain quantitative checks by applying the 
same interpretative scheme to the classical ether-drift experiments (Michelson-Morley, Miller, 
Illingworth, Joos,...). These old experiments were performed with interferometers where light 
was propagating in air or helium at atmospheric pressure. In this regime, where {M — 1) is 
a very small number, the theoretical fringe shifts expected on the basis of Eqs. (|23p and (|24p 
are much smaller than the classical prediction 0{/3'^) and it becomes conceivable that tiny 
non-zero effects might have been erroneously interpreted as 'null results'. 

To make this more evident, let us adopt Eq. (|24p . Then, an anisotropy of the two-way 
velocity of light could be measured by rotating a Michelson interferometer. As anticipated, in 
the rest frame of the apparatus, the length L of its arms does not depend on their orientation 
so that the interference pattern between two orthogonal beams of light depends on the time 
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difference 



AT{e) 



2L 2L 



(28) 



In this way, by introducing the wavelength A of tlie light source and the projection v of the 

2 

relative velocity in the plane of the interferometer, one finds to order ^ the fringe shift 



In the above equation the angle 9q = 9o{t) indicates the apparent direction of the ether-drift 
in the plane of the interferometer (the 'azimuth') and the square of the observable velocity 



is re-scaled by the very small quantity 2(J\f — 1) with respect to the true kinematical velocity 
f^(t). Notice that relation ()29p is formally identical to the classical prediction of a second- 
harmonic effect with the only replacement v — t- Vohs- 

Therefore, in this scheme, the interpretation of the experiments is transparent. According 
to Special Relativity, there can be no fringe shift upon rotation of the interferometer. In 
fact, if light propagates in a medium, the frame of isotropic propagation is always assumed 
to coincide with the laboratory frame, where the container of the medium is at rest, and thus 
one has fobs = v = 0. On the other hand, if there were fringe shifts, one could try to deduce 
the existence of a preferred frame S 7^ S" provided the following minimal requirements are 
fulfilled : i) the fringe shifts exhibit an angular dependence of the type in Eq. (j29p ii) by using 
gaseous media with different refractive index one gets consistency with Eq. (j30p in such a way 
that different fobs correspond to the same kinematical v. 

Before starting with the analysis of the classical experiments, one more remark is in order. 
In principle, even a single observation, within its experimental accuracy, can determine the 
existence of an ether-drift. However interpretative models are required to compare results 
obtained at different times and in different places. In the scheme of Eqs. (|29p and (j30p . the 
crucial information is contained in the two time-dependent functions v = v{t) and = ^o(*)) 
respectively the magnitude of the velocity and the apparent direction of the azimuth in the 
plane of the interferometer. For their determination, the standard assumption is to consider 
a cosmic Earth's velocity with well defined magnitude V , right ascension a and angular 
declination 7 that can be considered constant for short-time observations of a few days where 
there are no appreciable changes due to the Earth's orbital velocity around the Sun. In this 
framework, where the only time dependence is due to the Earth's rotation, one identifies 




(29) 




(30) 
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v{t) = v{t) and 6o{t) = OQ{t) where v{t) and 6o{t) derive from the simple application of 
spherical trigonometry [55] 

cos z{t) = sin 7 sin cp + cos 7 cos cj) cos(r — a) (31) 

^'•^^ ^ — sin z(t) cos 6o{t) = sin ^ cos (j) — cos J sin 4>cos{t — a) (32) 



V 

_ g-j^ ^^^-j g-j^ 0Q[t) = COS 7 sin(r — a) (33) 



^'(O = ^v^{t) + v^it) = V sin z{t), (34) 

Here z = z{t) is the zenithal distance of V, </> is the latitude of the observatory, r = Wgidi 
is the sidereal time of the observation in degrees (wsid ~ 2^&) ^"^^ angle 9q is counted 
conventionally from North through East so that North is 0o = and East is = 90°. 
To explore the observable implications, let us first re- write the basic Eq. ()29p as 

AA(0) 2L{N -l)v'^{t) 



X A c2 

where 



■ cos 2(9 - Ooit)) = 2C{t) cos 29 + 2S{t) sin 29 (35) 



L{N-l) v^{t) L{M-l) v^{t) . 

C(t) = ^ cos 26^0 (*) = 5— sm26'o(*) (36) 

A A 



Then Eqs. (|3T|) — p4|) amount to the structure 

S{t) = S{t) = Ssi sin T + Sci cos t + Ss2 sin(2r) + Sc2 cos(2r) (37) 

C{t) = C{t) = Co + C^i sin T + Cci cos r + C^2 sin(2r) + (7^2 cos(2t) (38) 
with Fourier coefficients iJZ = — ^■^) 

Co = -^7^(3cos27 - l)cos2(/> (39) 

Csi = — — 7^sinQsin27sin2(/) Cci = — — '7^cosasin27sin2(/» (40) 

C^2 = ^■7^sin2acos^7(l + sin^(/)) Cc2 = ^7e.cos 2q cos^ 7(1 + sin^ 0) (41) 

and 

5.1 = --^ Sci = ^ (42) 

sm (p sm <p 

_ 2sin,/. _ 2sin0 

1 + sm d) 1 + sm m 
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These standard forms are nowadays adopted in the analysis of the data of the ether-drift 
experiments j46j . However, one should not forget that Eq. (j24|) represents only an approxi- 
mation for the full structure Eq. (j23p . Therefore, even for short-time observations, one might 
not obtain from the data completely consistent determinations of the kinematical parameters 
{V, a, 7). In addition, by using a physical analogy, and by representing the physical vacuum as 
a fluid, the scheme Eqs. (|37p . (|38p of smooth sinusoidal variations associated with the Earth's 
rotation corresponds to the conditions of a pure laminar flow associated with simple regular 
motions. Instead, the physical vacuum might be similar to a form of turbulent fluid, where 
large-scale and small-scale flows are only indirectly related. 

In this modified perspective, which finds motivations in some basic foundational aspects 
of both quantum physics and relativity [56l EZl [581 EHl EO] and in those representations of the 
vacuum as a form of 'space-time foam' which indeed resembles a turbulent fiuid |14 t ll5 [[TUlll7| . 
the ether-drift might exhibit forms of time modulations that do not fit in the scheme of 
Eqs. (j37p .(|38 p . To evaluate the potential effects, and by still retaining the functional form 
Eq. ([5?]) . one could first re- write Eqs.([3HD as 

C(t) = ^(-^"^^ ^'W-^'(^) ^(^) ^ MAA-l) 2v,it)vy{t) 

A A 

where Vx{t) = v(t) cos6o(t) and Vy{t) = v{t) sinOo^t). Then, by exploiting the turbulence 

scenario, one could model the two velocity components Vx{t) and Vy(t) as stochastic fiuctua- 

tions. In this different scheme, where now v{t) ^ v{t) and 9o{t) 7^ Oo{t), experimental results 

which, on consecutive days and at the same sidereal time, deviate from Eqs. (|3ip — (j34p do not 

necessarily represent spurious effects. Equivalently, if data collected at the same sidereal time 

average to zero this does not necessarily mean that there is no ether-drift. This particular 

aspect will be discussed at length in the rest of the paper. 

After this important premise, we shall now proceed in Sects. 3-8 with our re-analysis of 

the classical experiments. In the end. Sect. 9 will contain a summary, a brief discussion of the 

modern experiments and our conclusions. 



3. The original Michelson-Morley experiment 

The Michelson-Morley experiment |61j is probably the most celebrated experiment in the 
history of physics. Its result and its interpretation have been (and are still) the subject of 
endless controversies. For instance, for some time there was the idea [62] that, by taking into 
account the reflection from a moving mirror and other effects, the predicted shifts would be 
largely reduced and become unobservable. These points of view are summarized in Hedrick's 
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contribution to the 'Conference on the Michelson-Morley experiment' [63] (Pasadena, Febru- 
ary 1927) which was attended by the greatest experts of the time, in particular Lorentz and 
Michelson. The arguments presented by Hedrick were, however, refuted by Kennedy [64J in a 
paper of 1935 where, by using Huygens principle, he re-obtained to order jt? the classical 
result of Eq. ()29p (with the identification fobs = ^')- 




Figure 1: The Michelson-Morley fringe shifts as reported by Hicks [6^1 . Solid and dashed lines 
refer respectively to noon and evening observations. 

In this framework, the fringe shift is a second-harmonic effect, i.e. periodic in the range 
[0, vr], whose amplitude A2 is predicted differently by using the classical formulas or Lorentz 
transformations ()29p 

^class ^LtP 1 ^ LVoh^ ^ _ ^x^class ^^^^ 

Notice also that upon rotation of 7r/2 with respect to 9 = the predicted fringe shift is 2^2- 
Now, for the Michelson-Morley interferometer the whole effective optical path was about 



14 



L = 11 meters, or about 2 • 10 in units of light wavelengths, so for a velocity f ~ 30 km/s 
(the Earth's orbital velocity about the Sun, and consequently the minimum anticipated drift 
velocity) the expected classical 2nd-harmonic amplitude was A^i'^"" ~ 0.2. This value can thus 
be used as a reference point to obtain an observable velocity, in the plane of the interferometer, 
from the actual measured value of A2 through the relation 



Michelson and Morley performed their six observations in 1887, on July 8th, 9th, 11th and 
12th, at noon and in the evening, in the basement of the Case Western University of Cleveland. 
Each experimental session consisted of six turns of the interferometer performed in about 36 
minutes. As well summarized by Miller in 1933 [HS], "The brief series of observations was 
sufficient to show clearly that the effect did not have the anticipated magnitude. However, 
and this fact must be emphasized, the indicated effect was not zero" . 

The same conclusion had already been obtained by Hicks in 1902 [66]: " ..the data published 
by Michelson and Morley, instead of giving a null result, show distinct evidence for an effect 
of the kind to be expected". Namely, there was a second-harmonic effect. But its amplitude 
was substantially smaller than the classical expectation (see Fig.l). 

Quantitatively, the situation can be summarized in Figure 2, taken from Miller [65], where 
the values of the effective velocity measured in various ether-drift experiments are reported 
and compared with a smooth curve fitted by Miller to his own results as function of the 
sidereal time. 

For the Michelson-Morley experiment, the average observable velocity reported by Miller 
is about 8.4 km/s. Comparing with the classical prediction for a velocity of 30 km/s, this 
means an experimental 2nd- harmonic amplitude 



which is about twelve times smaller than the expected result. 

Neither Hicks nor Miller reported an estimate of the error on the 2nd harmonic extracted 
from the Michelson-Morley data. To understand the precision of their readings, we can look 
at the original paper [61] where one finds the following statement : "The readings are divisions 
of the screw- heads. The width of the fringes varied from 40 to 60 divisions, the mean value 
being near 50, so that one division means 0.02 wavelength". Now, in their tables Michelson 
and Morley reported the readings with an accuracy of 1/10 of a division (example 44.7, 44.0, 
43.5,..). This means that the nominal accuracy of the readings was ±0.002 wavelengths. In 




(46) 




0.2 (_)2^ 0.016 



(47) 
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Fig. 4. Velocity of ether drift observed by Michelson 
and Morley in and by Morby and Miller in 1902, 

lfK)4 and 1905, compared with the velot^ity obtained by 
Miller in 1925. 



Figure 2: The magnitude of the observable velocity measured in various experiments as reported 
by Miller 165^ . 

fact, in units of wavelengths, they reported values such as 0.862, 0.832, 0.824,.. Furthermore, 
this estimate of the error agrees well with Born's book [67J. In fact. Born, when discussing 
the classically expected fractional fringe shift upon rotation of the apparatus by 90°, about 
0.37, explicitly says: "Michelson was certain that the one-hundredth part of this displacement 
would still be observable" (i.e. 0.0037). Therefore, to be consistent with both the original 
Michelson-Morley article and Born's quotation of Michelson's thought, we shall adopt ±0.004 
as an estimate of the error 0. 

With this premise, the Michelson-Morley data were re-analyzed in ref.|51]. To this end, 
one should first follow the well defined procedure adopted in the classical experiments as 
described in Miller's paper [65]. Namely, by starting from each set of seventeen entries (one 
every 22.5°), say E{i), one has first to correct the data for the observed linear thermal drift. 
This is responsible for the difference E{1) — E{17) between the 1st entry and the 17th entry 
obtained after a complete rotation of the apparatus. In this way, by adding 15/16 of the 
correction to the 16th entry, 14/16 to the 15th entry and so on, one obtains a set of 16 

''To confirm that such estimate should not be considered unreafistically small, we report explicitly Michel- 
son's words from ref. [6^ : "I must say that every beginner thinks himself lucky if he is able to observe a shift 
of 1/20 of a fringe. It should be mentioned however that with some practice shifts of 1/100 of a fringe can be 
measured, and that in very favorable cases even a shift of 1/1000 of a fringe may be observed." 
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Table 1: The fringe shifts — for all noon (n.) and evening (e.) sessions of the Michelson- 
Morley experiment. 



i 


July 8 (n.) 


July 9 (n.) 


July 11 (n.) 


July 8 (e.) 


July 9 (e.) 


July 12 (e.) 


1 


-0.001 


+0.018 


+0.016 


-0.016 


+0.007 


+0.036 


2 


+0.024 


-0.004 


-0.034 


+0.008 


-0.015 


+0.044 


3 


+0.053 


-0.004 


-0.038 


-0.010 


+0.006 


+0.047 


4 


+0.015 


-0.003 


-0.066 


+0.070 


+0.004 


+0.027 


5 


-0.036 


-0.031 


-0.042 


+0.041 


+0.027 


-0.002 


6 


-0.007 


-0.020 


-0.014 


+0.055 


+0.015 


-0.012 


7 


+0.024 


-0.025 


+0.000 


+0.057 


-0.022 


+0.007 


8 


+0.026 


-0.021 


+0.028 


+0.029 


-0.036 


-0.011 


9 


-0.021 


-0.049 


+0.002 


-0.005 


-0.033 


-0.028 


10 


-0.022 


-0.032 


-0.010 


+0.023 


+0.001 


-0.064 


11 


-0.031 


+0.001 


-0.004 


+0.005 


-0.008 


-0.091 


12 


-0.005 


+0.012 


+0.012 


-0.030 


-0.014 


-0.057 


13 


-0.024 


+0.041 


+0.048 


-0.034 


-0.007 


-0.038 


14 


-0.017 


+0.042 


+0.054 


-0.052 


+0.015 


+0.040 


15 


-0.002 


+0.070 


+0.038 


-0.084 


+0.026 


+0.059 


16 


+0.022 


-0.005 


+0.006 


-0.062 


+0.024 


+0.043 



corrected entries 

Ecorrii) = ^(^(1) - ^(17)) + E{i) (48) 
16 

The fringe shifts are then defined by the differences between each of the corrected entries 
Ecovvii) and their average value (-Ecorr) as 

= ^corr(i) - (^corr) (49) 

The resulting data are reported in Table 1. 

With this procedure, the fringe shifts Eq. ()49p are given as a periodic function, with van- 
ishing mean, in the range < < 27r, with 9 = 27r, so that they can be reproduced in a 
Fourier expansion. Notice that in the evening observations the apparatus was rotated in the 
opposite direction to that of noon. 
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Table 2: The amplitude of the fitted second- harmonic component ^f'^^ for the six experimental 
sessions of the Michelson-Morley experiment. 



SESSION 
July 8 (noon) 
July 9 (noon) 
July 11 (noon) 
July 8 (evening) 
July 9 (evening) 
July 12 (evening) 



0.010 ± 0.005 
0.015 ± 0.005 
0.025 ± 0.005 
0.014 ± 0.005 
0.011 ±0.005 
0.024 ± 0.005 



One can thus extract the amplitude and the phase of the 2nd-harmonic component by 
fitting the even combination of fringe shifts 

2A 

(see Fig. 3). This is essential to cancel the Ist-harmonic contribution originally pointed out 
by Hicks [66]. Its theoretical interpretation is in terms of the arrangements of the mirrors 
and, as such, this effect has to show up in the outcome of real experiments. For more details, 
see the discussion given by Miller, in particular Fig. 30 of ref.[65j, where it is shown that his 
observations were well consistent with Hicks' theoretical study. The observed Ist-harmonic 
effect is sizeable, of comparable magnitude or even larger than the second-harmonic effect. 
The same conclusion was also obtained by Shankland et al. |68j in their re-analysis of Miller's 
data. The 2nd-harmonic amplitudes from the six individual sessions are reported in Table 2. 

Due to their reasonable statistical consistency, one can compute the mean and variance of 
the six determinations reported in Table 2 by obtaining A^^^ ~ 0.016 ± 0.006. This value is 
consistent with an observable velocity 

?;obs ~ 8.4l|:^ km/s (51) 

Then, by using Eq. (j30|) . which connects the observable velocity to the projection of the kine- 
matical velocity in the plane of the interferometer through the refractive index of the medium 
where light propagation takes place (in our case air where N ~ 1.00029), we can deduce the 
average value 

V ~ 3491^^ km/s (52) 
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July 1 1 noon 

0.04 1 , 






Figure 3: A fit to the even combination B{9) Eg. /i50\) . The second harmonic amplitude is 
^EXP ^ 0.025 ± 0.005 and the fourth harmonic is A^^^ = 0.004 ± 0.005. The figure is taken 
from ref.^SJJ. Compare the data with the solid curve of July 11th shown in Fig.l. 

While the individual values of A2 show a reasonable consistency, there are substantial 
changes in the apparent direction of the ether-drift effect in the plane of the interferometer. 
This is the reason for the strong cancelations obtained when fitting together all noon sessions 
or all evening sessions [69]. For instance, for the noon sessions, by taking into account that the 
azimuth is always defined up to ±180°, one choice for the experimental azimuths is 357° it 14°, 
285° ± 10° and 317° ± 8° respectively for July 8th, 9th and 11th. For this assignment, the 
individual velocity vectors Vobs(cos Oq, — sin ^o) and their mean are shown in Fig. 4. According 
to the usual interpretation, the large spread of the azimuths is taken as indication that any 
non-zero fringe shift is due to pure instrumental effects. However, as anticipated in Sect. 2, 
this type of discrepancy could also indicate an unconventional form of ether-drift where there 
are substantial deviations from Eq. (j24p and/or from the smooth trend in Eqs. (j3ip — (j34p . For 
instance, in agreement with the general structure Eq. (|23p . and differently from July 11 noon, 
which represents a very clean indication, there are sizeable 4th- harmonic contributions (here 
^EXP ^ Q_oi9 ± 0.005 and A^^^ = 0.008 ± 0.005 for the noon sessions of July 8 and July 9 
respectively). In any case, the observed strong variations of are in qualitative agreement 
with the analogous values reported by Miller. To this end, compare with Fig. 22 of ref.[65] and 
in particular with the large scatter of the data taken around August 1st, as this represents the 
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Figure 4: The observable velocities for the three noon sessions and their mean. The x-axis 
corresponds to 9^ = = 360° and the y-axis to 9q = 270°. Statistical uncertainties of the 
various determinations are ignored. All individual directions could also be reversed by ISOP . 

epoch of the year which is closer to the period of July when the Michelson-Morley observations 
were actually performed. Thus one could conclude that individual experimental sessions 
indicate a definite non-zero ether-drift but the azimuth does not exhibit the smooth trend 
expected from the conventional picture Eqs. (j3T|) — (plj) . 

For completeness, we add that the large spread of the 0o~values might also reflect a 
particular systematic effect pointed out by Hicks |66j . As described by Miller [65], " before 
beginning observations the end mirror on the telescope arm is very carefully adjusted to secure 
vertical fringes of suitable width. There are two adjustments of the angle of this mirror which 
will give fringes of the same width but which produce opposite displacements of the fringes 
for the same change in one of the light-paths" . Since the relevant shifts are extremely small, 
"...the adjustments of the mirrors can easily change from one type to the other on consecutive 
days. It follows that averaging the results of different days in the usual manner is not allowable 
unless the types are all the same. If this is not attended to, the average displacement may 
be expected to come out zero — at least if a large number are averaged" [66j. Therefore 
averaging the fringe shifts from various sessions represents a delicate issue and can introduce 
uncontrolled errors. Clearly, this relative sign does not affect the values of A2 and this is 
why averaging the 2nd-harmonic amplitudes is a safer procedure. However, it can introduce 
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spurious changes in the apparent direction 9q of the ether-drift. In fact, an overall change of 
sign of the fringe shifts at all ^—values is equivalent to replacing ~^ i 7r/2. As a matter 
of fact, Hicks concluded that the fringes of July 8th were of different type from those of the 
remaining days. Thus for his averages (in our Fig.l) "the values of the ordinates are one-third 
of July 9 + July 11 — July 8 and one-third of July 9 + July 12 — July 8" f66] for noon and 
evening sessions respectively. If this were true, one choice for the azimuth of July 8th could 
now be e^^^ = 267° ± 14°. This would orient the arrow of July 8th in Fig. 4 in the direction 
of the y— axis and change the average azimuth from {9q^^) ~ 317° to {9^'^^) ~ 290°. We'll 
return to this particular aspect in our Appendix II. 

Let us finally compare with the interpretation that Michelson and Morley gave of their 
data. They start from the observation that " ...the displacement to be expected was 0.4 fringe" 
while "...the actual displacement was certainly less than the twentieth part of this". In this 
way, since the displacement is proportional to the square of the velocity, "...the relative velocity 
of the earth and the ether is... certainly less than one-fourth of the orbital earth's velocity". 
The straightforward translation of this upper bound is Vobs < 7.5 km/s. However, this estimate 
is likely affected by a theoretical uncertainty. In fact, in their Fig. 6, Michelson and Morley 
reported their measured fringe shifts together with the plot of a theoretical second-harmonic 
component. In doing so, they plotted a wave of amplitude A2 = 0.05, that they interpret 
as one-eight of the theoretical displacement expected on the base of classical physics, thus 
implicitly assuming A2^^^=0A. As discussed above, the amplitude of the classically expected 
second-harmonic component is not 0.4 but is just one-half of that, i.e. 0.2. Therefore, their 
experimental upper bound A'^^ < ^ =0.02 is actually equivalent to fobs < 9.5 km/s. If 
we now consider that their estimates were obtained after superimposing the fringe shifts 
obtained from various sessions (where the overall effect is reduced, see our Fig.l), we deduce 
a substantial agreement with our result Eq. ([?T|) . 

4. Morley-Miller 

After the original 1887 experiment, there was much interest in the Michelson-Morley result 
that, being too small to meet any classical prediction, was apparently contradicting two 
cornerstones of physics: Galilei's transformations and/or the existence of the ether. For this 
reason, one of the most influential physicists of the time. Lord Kelvin, after his conference at 
the 1900 Paris Expo, induced Morley and his young collaborator Dayton Miller to design a 
new interferometer (where the effective optical path was increased up to 32 meters) to improve 
the accuracy of the measurement over the 1887 result. 
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It must be emphasized that Morley and Miher [70], in their observations of 1905, super- 
imposed the data of the morning with those of the evening. As explained by Miher [B^ , the 
two physicists were assuming that the ether drift had to be obtained by combining the motion 
of the solar system relative to nearby stars, i.e. toward the constellation of Hercules with a 
velocity of about 19 km/s, with the annual orbital motion ("We now computed the direction 
and the velocity of the motion of the centre of the apparatus by compounding the annual 
motion in the orbit of the earth with the motion of the solar system toward a certain point 
in the heavens... There are two hours in each day when the motion is in the desired plane 
of the interferometer" [70]). The observations at the two times (about 11:30 a.m. and 9:00 
p.m.) were, therefore, combined in such a way that the presumed azimuth for the morning 
observations coincided with that for the evening ( "The direction of the motion with reference 
to a fixed line on the floor of the room being computed for the two hours, we were able to 
superimpose those observations which coincided with the line of drift for the two hours of 
observation" [70]). However, the observations for the two times of the day gave results having 
nearly opposite phases. When these were combined, the result was nearly zero. For this 
reason, the value then reported of an observable velocity of 3.5 km/s is incorrect and does 
not correspond to the actual results of the basic observations. The error was later understood 
and corrected by Miller who found that the two sets of data were each indicating an effective 
velocity of about 7.5 km/s (see Figure 11 of Miller's paper |65j). For this reason, the correct 
average observable velocities for the entire period 1902-1905 are those shown in our Figure 2 
between 7 and 10 km/s or 

Vohs ~ (8.5 ± 1.5) km/s (53) 
By using Eq. pO]) . we then deduce the average value 

V ~ (353 ± 62) km/s (54) 

5. Kennedy-Illingworth 

An interesting development was proposed by Kennedy in 1926. As summarized in his con- 
tribution to the previously mentioned Conference on the Michelson- Morley experiment [63], 
his small optical system was enclosed in an effectively insulated, sealed metal case containing 
helium at atmospheric pressure. Because of its small size, "...circulation and variation in 
density of the gas in the light paths were nearly eliminated. Furthermore, since the value of 
M — 1 is only about 1/10 that for the air at the same pressure, the disturbing changes in 
density of the gas correspond to those in air to only 1/10 of the atmospheric pressure". The 
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Table 3: The infra-session averages {Da) and (Db) obtained from the 10 sets of rotations in 
each of the 32 sessions of lUingworth's experiment. These values have been obtained from the 
weights of lUingworth's Table III by applying the conversion factor 0.002. 



5 A.M. 


5 A.M. 


11 A.M. 


11 A.M. 


5 P.M. 


5 P.M. 


11 P.M. 


11 P.M. 


(Da) 


(Db) 


(Da) 


(Db) 


(Da) 


(Db) 


(Da) 


(Db) 


+0.00024 


-0.00066 


+0.00070 


-0.00022 


+0.00024 


+0.00044 


-0.00010 


+0.00024 


+0.00114 


+0.00024 


-0.00042 


-0.00036 


-0.00056 


-0.00046 


+0.00018 


+0.00018 


+0.00000 


+0.00000 


-0.00006 


-0.00052 


-0.00144 


-0.00080 


-0.00126 


-0.00006 


+0.00020 


-0.00044 


-0.00030 


+0.00012 


-0.00016 


+0.00004 


-0.00044 


-0.00026 


+0.00064 


+0.00000 


-0.00022 


+0.00038 


+0.00018 


+0.00016 


+0.00000 


+0.00024 


-0.00002 


-0.00010 


+0.00048 


+0.00020 


+0.00030 


+0.00030 


-0.00040 


-0.00004 






-0.00014 


-0.00006 


+0.00030 


+0.00014 










-0.00006 


+0.00004 


+0.00036 


-0.00036 










-0.00006 


+0.00016 


+0.00006 


-0.00006 










+0.00000 


+0.00024 


-0.00010 


+0.00010 







essential ingredient of Kennedy's apparatus consisted in the introduction of a small step, 1 /20 
of wavelength thick, in one of the total reflecting mirrors of the interferometer allowing, in 
principle, for an ultimate fringe shift accuracy 1 • 10~^. To take full advantage of this possibil- 
ity, Kennedy should have disposed of perfect mirrors and of a suitable (hotter) source of light. 
In the original version of the experiment, these refinements were not implemented giving an 
actual fringe shift accuracy of 2 • 10"^. In these conditions, as Kennedy explicitly savs|63j. 
"...the velocity of 10 km/s found by Prof. Miller would produce a fringe shift corresponding 
to 8 • 10"'^", four times larger than the experimental resolution. Since the effect is quadratic 
in the velocity, Kennedy's result, fringe shifts < 2 • 10"^, can then be summarized as 

fobs < 5 km/s (55) 

By using Eq. (j30p . for helium at atmospheric pressure where M ~ 1.000035, this bound 
amounts to restrict the kinematical value hy v < 600 km/s. 

Kennedy's apparatus was further refined by Illingworth in 1927 [7T]. Besides improving 
the quality of the mirrors and of the source, lUingworth's data taking was also designed to 
reduce the presence of steady thermal drift and of odd harmonics. Looking at lUingworth's 
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paper, one finds that his refinements reached indeed the nominal 0(10"^) accuracy mentioned 
by Kennedy, namely about 1/1500 of wavelength for the individual readings and (1-^2)- 10"'^ 
at the level of average values. 

Let us now analyze Illingworth's results. He performed four series of observations in the 
first ten days of July 1927. These consisted of 32 experimental sessions, conducted daily 
at 5 A.M. (6), 11 A.M. (10), 5 P.M. (10) and 11 P.M.(6), in which he was measuring the 
fringe displacement caused by a rotation through a right angle of the apparatus. To take into 
account 90° rotations let us first re- write Eq. ()29p as 

^^=A2Cos2{e-eo) (56) 
A 

Therefore Illingworth, in his first set (set A) of 10 rotations. North, East, South, West and 
back to North, was actually measuring Da = 2yl2 cos 26o. In a second set (set B), North-East, 
North- West, South- West, South-East and back to North-East, performed immediately after 
the set A, he was then measuring Db = 2A2 sin 2^0- Notice that both Da and Db differ from 
the positive-definite quantity D = 2A2 that should be inserted in Illingworth's numerical 
relation for his apparatus Vohs = 112-v/l). Therefore, the reported values for the two velocities 
VA = 112-y/|L)A| and vb = 112a/|I?_b| should only be taken as lower bounds for the true Vohs- 
The mean values {Da) and (Db) obtained from the 10 sets of rotations in the 32 individual 
sessions can be obtained from Illingworth's Table III and, for the convenience of the reader, 
are reported in our Table 3. 

From Table 3, one finds that the quantity {D a)"^ + {Db)'^ has a mean value of about 
0.00045, which corresponds to Vohs ~ 2.4 km/s. Thus, by using Eq. (j30p for helium at atmo- 
spheric pressure, we would tentatively deduce an average value v ~ 284 km/s. 

However, this is only a very partial view. To go deeper into Illingworth's experiment we 
have to consider his basic measurements, i.e. the individual turns of his interferometer. In 
this case, the only known basic set of data reported by Illingworth is set A of July 9th, 11 
A.M. This set has been re-analyzed by Munera [72j and his values for the fringe shifts are 
reported in our Table 4. 

As one can see, the fringe shifts are not small and correspond to an observable velocity in 
the range 2-5 km/s. However, their sign seems to change randomly. Therefore, if one attempts 
to extract the observable velocity from the mean of the 10 determinations, {Da) ~ —0.00006, 
the resulting value 0.9 km/s is much smaller than all individual determinations. The basis 
of Munera's analysis was instead to estimate fobs from from which he obtained an 

average velocity fobs = 3.13 it 1.04 km/s. 
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Table 4: Illingworth's set A of July 9th, 11 A.M. as re-analyzed by Munera J7^y . 



Rotation 


Da 


|DaI 


VA[km/s] 


1 


-0.00100 


+0.00100 


3.54 


2 


+0.00066 


+0.00066 


2.89 


3 


-0.00066 


+0.00066 


2.89 


4 


-0.00066 


+0.00066 


2.89 


5 


-0.00166 


+0.00166 


4.57 


6 


+0.00234 


+0.00234 


5.41 


7 


+0.00100 


+0.00100 


3.54 


8 


+0.00034 


+0.00034 


2.04 


9 


+0.00000 


+0.00000 


0.00 


10 


-0.00100 


+0.00100 


3.54 



Now, the standard interpretation of such apparently random changes of sign is in terms 
of typical instrumental effects and the standard method for eliminating these is the original 
averaging procedure as employed by Illingworth. But we will now show that they could also 
indicate an unconventional form of stochastic drift, of the type already mentioned in the 
previous sections, and in which Miinera's re-estimate has a definite significance. To this end, 
we shall first use the relations 

DA{t)=4Cit) DBit)=4S{t) (57) 

where the two functions C(t) and S(t) have been introduced in Eqs. ([36|) and ([H]) . Thus 
Eqs. (|57|) can be re- written as 

8L{M-1) vUt)-vlit) SLjM-l) v.jt) vyjt) 
L>A{t) - UB[t) - -2 (58) 



where Vx{t) = v{t) cosOoit) and Vy{t) = v{t)s\n.9Q{t). In this way, by using the numerical 

L 
A 



relation for Illingworth's experiment k- (^Q^™/'^) ^ 0.035 and the value of the helium refractive 



index, we obtain 

The required random ingredient can then be introduced by characterizing the two velocity 
components Vx{t) and Vy{t) as turbulent fluctuations. To this end, there can be several ways. 
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Here we shall restrict to the simplest choice of a turbulence which, at small scales, appears 
statistically isotropic and homogeneous @. This represents a zeroth-order approximation which 
is useful to illustrate basic phenomenological features associated with an underlying stochastic 
vacuum. To explore the resulting temporal pattern of the data, we have followed refs.|741 175j 
where velocity flows, in statistically isotropic and homogeneous 3-dimensional turbulence, are 
generated by unsteady random Fourier series. The perspective is that of an observer moving 
in the turbulent fluid who wants to simulate the two components of the velocity in his x-y 
plane at a given fixed location in his laboratory. This leads to the general expressions 

oo 

fJxit) = ^ COS u)nt + x„(2) sina;„t] (60) 

71=1 

oo 

= X] coswnt + y„(2) sin (61) 

n=l 

where a;„ = Inir/T, T being a time scale which represents a common period of all stochastic 
components. We have adopted the typical value T = Tday= 24 hours. However, we have 
also checked with a few runs that the statistical distributions of the various quantities do not 
change substantially by varying T in the rather wide range 0.1 Tday < T < 10 Tjay 

The coefficients Xn{i = 1,2) and y„(i = 1,2) are random variables with zero mean. They 
have the physical dimension of a velocity and we shall denote by [—v, v] the common interval 
for these four parameters. In terms of v the statistical average of the quadratic values can be 
expressed as 

{xlii = l,2))stat = {ylii = l,2))stat = ^ (62) 

for the uniform probability model (within the interval [— 'U,^]) which we have chosen for 
our simulations. Finally, the exponent ij controls the power spectrum of the fluctuating 
components. For the simulations, between the two values rj = 5/6 and rj = 1 reported in 
ref.[75], we have chosen rj = 1 which corresponds to the point of view of an observer moving 
in the fluid. 

Thus, within this simple model for D^it) and Dsit), v is the only parameter whose 
numerical value could reflect the properties of a large-scale motion, for instance of the Earth's 
motion with respect to the Cosmic Microwave Background (CMB). For this reason, here, we 
have adopted the fixed value v = Vcmb = 370 km/s. With these premises, our results can be 
illustrated by first considering the basic set of 10 complete rotations of the apparatus during 
which Illingworth's fringe shifts (produced by 90° rotations) were recorded every 30 seconds. 

'This picture reflects the basic Kolmogorov theory .73.^ of a fluid with vanishingly smaU viscosity. 
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Figure 5: A simulation o/DA(t) and DB(t), in units 10"^ and every 30 seconds, from typical 
sequences of 1200 seconds. The average values are {Da) = 0.00028 and (Db) = 0.00011. The 
velocity parameter is v = Vcmb = 370 km/s. 

Therefore, this type of simulations consists in generating 40 values during a total time of 1200 
seconds. As an illustration, two typical sequences of DA(t) and DB(t), in units 10"'^, are 
shown in Fig. 5. 

As one can see, the magnitude 0(10"^) and the random nature of the instantaneous values 
is completely consistent with the entries of Table 4. Also the resulting infra-session averages 
(Da) = 0.00028 and (Db) = 0.00011 are completely consistent with the typical entries of 
Table 3. 

To obtain further insight, we have then performed extensive simulations for large sequences 
of measurements. The histograms of a set of 10000 determinations of -D^(t) and DB{t) (again 
generated every 30 seconds) are reported in panels (a) and (b) of Fig. 6. 

Notice that these distributions are clearly "fat-tailed" and very different from a Gaussian 
shape. This kind of behavior is characteristic of probability distributions for instantaneous 
data in turbulent flows (see e.g. |761 177j). To better appreciate the deviation from Gaussian 
behavior, in panels (c) and (d) we plot the same data in a log— log scale. The resulting 
distributions are well fitted by the so-called g— exponential function |7B] 



with entropic index q ~ 1.1. For such large samples of data, the statistical averages {Da) and 



f,{x)=a{l-{l-q)xbf'^^-'i^ 



(63) 
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Figure 6: We show, see (a) and (h), the histograms W obtained from a simulation for 
Da = DA(t) and Db = DB(t). The vertical normalization is to a unit area. The mean values 
are (Da) = 0.75 • 10~^, (Db) = -1.1 • 10~^ and the standard deviations ct(Da) = 0.75 • 10~^, 
(t(Db) = 0.83 • 10^'^. We also show, see (c) and (d), the corresponding plots in logarithmic 
scale and the fits with Eg. h63^) . The parameters of the fit are q=1.07, a=2 and h=2.2 for Da 
and q=1.12, a=2 and h=2.3 forY)-Q. The total statistics correspond to 10.000 values generated 
at steps of 30 seconds. The velocity parameter is v = Vcmb = 370 km/s. 
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{Db) are vanishingly small in units of the typical instantaneous values 0(10"^) and any non- 
zero average has to be considered as statistical fluctuation. On the other hand, the standard 
deviations a{DA) and a{DB) have definite non-zero values which reflect the magnitude of the 
scale parameter v. By keeping v fixed at 370 km/s, we have found 

aiDA) ~ (0.74 ± 0.05) • 10"^ (t{Db) ~ (0.83 ± 0.06) • 10"^ (64) 

whose uncertainties reflect the observed variations due to the truncation of the Fourier modes 
in Eqs. (f60|) . ([6T]) and to the dependence on the random sequence. Taking this calculation 
into account gives a mean spread slightly less, about 0.65 • 10~^, for the effect of stochastic 
drift in Illingworth's measurements. This is comparable to the uncertainty of the individual 
readings which, in the best case, was of 1/1500 wavelengths, i.e. ±0.7 • 10^^. By combining 
in quadrature the two uncertainties, one gets a good agreement with our Table 4 where the 
variance of the mean is about ±1 • 10~^. Finally, the simulation is also useful to get indications 
on the expected value of the observable velocity. In fact, with vanishingly small values of {Da_) 
and (Db) one gets (-D^) ~ cr^(D^) and {D'^) ~ g'^^Db)- Therefore one obtains the following 
two average estimates of Vohs 

Vohs ~ ll2^/a{DX) ~ 3.05 km/s v^hs ~ ll2^/a{DB) ~ 3.23 km/s (65) 

with a mean value of 3.14 km/s which is very close to Munera's determination fobs = 3.13ibl.04 
km/s. 

We emphasize that one could further improve the stochastic model by introducing time 
modulations and/or slight deviations from isotropy. For instance, v could become a function 
of time V = v[t). By still retaining statistical isotropy, this could be used to simulate the 
possible modulations of the projection of the Earth's velocity in the plane of the interferometer. 
Or, one could fix a range, say [— for the two random parameters 2;„(1) and Xn(2), 
which is different from the range [—Vy^Vy] for the other two parameters yn(l) and yn{2). 
Finally, Vx and Vy could also become given functions of time, for instance Vx{t) = v{t) cos ^o(^) 
Vy{t) = ■u(t) sin6'o(t), v{t) and 6o{t) being defined in Eqs. ([HT]) — P^ . We shall discuss this 
other alternative later on, in connection with Joos' 1930 experiment. 

In any case, by accepting this type of picture of the ether-drift, it is clear that further 
reduction of the data by performing inter-session averages (((...))) among the various sessions, 
can wash out completely the physical information contained in the original observations. In 
Table 5, we report the final inter-session averages {{Da_)) and {{Db)) obtained by Illingworth 
for the various observation times. 
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Table 5: lUingworth's final inter-session averages. 



Observations 


{{Da)) 


{{Db)) 


5 A.M. 


+0.00036 ± 0.00012 


-0.00016 ± 0.00009 


11 A.M. 


-0.00001 ± 0.00007 


-0.00000 ± 0.00006 


5 P.M. 


-0.00008 ± 0.00012 


-0.00005 ± 0.00008 


11 P.M. 


-0.00034 ± 0.00014 


+0.00005 ± 0.00006 



Nevertheless, in spite of the strong cancelations expected from the averaging reduction 
process mentioned above, some non-zero value is still surviving. Therefore, regardless of our 
simulations, one could draw the following conclusions. Traditionally, from these final averages 
for {{Da)) at 5 A.M. and at 11 P.M. one has been deducing the values va ~ 2.12 km/s and 
VA ~ 2.07 km/s respectively. Therefore, from these two estimates of va that, as anticipated, 
represent lower bounds for Wobsj it follows that there were values of fobs which clearly had to 
be larger than both. For this reason, it is surprising that in the literature (see for instance [68]) 
this 2.1 km/s velocity value reported by Illingworth, rather than being interpreted correctly as 
a /ow;er bound, was always interpreted incorrectly as an upper bound placed by his experiment. 
Instead, by combining with the previous Kennedy's upper bound t7obs < 5 km/s, one deduces 
that these two experiments, where light was propagating in helium at atmospheric pressure, 
give a range for the observable velocity 

(Kennedy + Illingworth) 2 km/s < t^obs < 5 km/s (66) 

in complete agreement with Miinera's determination 

Vohs = 3.1 + 1.0 km/s (67) 

Prom this last estimate, by using Eq. ()30p and taking into account that for helium at atmo- 
spheric pressure the refractive index is A/" ~ 1.000035, one obtains a kinematical velocity 

V ~ (370 + 120) km/s (68) 

consistently with the velocity values Eqs. ()52p and ()54p from the Michelson-Morley and Mor ley- 
Miller experiments. 
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6. Miller 



Munera's analysis j72| is also interesting because he applied the same method used for Illing- 
worth's observations to the only known Miller set of data explicitly reported in the literature. 
In this case, his value Vobs = 8.2 it 1.4 km/s, after correcting with Eq. (j30p . confirms the 
estimate v ~ 350 km/s for the average velocity in the plane of the interferometer. 

This close agreement with the Michelson-Morley value 8.4 km/s is also confirmed by the 
critical re-analysis of Shankland et al. [BH]. Differently from the original Michelson-Morley 
experiment, Miller's data were taken over the entire day and in four epochs of the year. 
However, after the critical re-analysis of the original raw data performed by the Shankland 
team, there is now an independent estimate of the average determinations A^^^ for the four 
epochs. Their values 0.042, 0.049, 0.038 and 0.045, respectively for April 1925, July 1925, 
September 1925 and February 1926 (see page 170 of ref.[UH]) are so well statistically consistent 
that one can easily average them. The overall determination from Table III of |68j 

= 0.044 ± 0.022 (69) 

when compared with the equivalent classical prediction for Miller's interferometer ^2^^'^'' = 
L (30km/s) ^ g corresponds to an average observable velocity 

Vohs = 8.4+^1 km/s (70) 

and, by using Eq. ([30]) . to a true kinematical value 

V = 349+^^4 km/s (71) 

We are aware that our conclusion goes against the widely spread belief, originating precisely 
from the paper of Shankland et al. ref . |68j . that Miller's results might actually have been due 
to statistical fluctuation and/or local temperature conditions. To a closer look, however, the 
arguments of Shankland et al. are not so solid as they appear when reading the Abstract of 
their paper [^. In fact, within the paper these authors say that "...there can be little doubt 
that statistical fluctuations alone cannot account for the periodic fringe shifts observed by 
Miller" (see page 171 of ref.|68)). Further, although "...there is obviously considerable scatter 
in the data at each azimuth position,... the average values. ..show a marked second harmonic 
effect" (see page 171 of ref. [68]). In any case, interpreting the observed effects on the basis 
of the local temperature conditions is certainly not the only explanation since "...we must 
admit that a direct and general quantitative correlation between amplitude and phase of the 



detailed rebuttal of the criticism raised by the Shankland team can be found in ref. [79) 
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observed second harmonic on the one hand and the thermal conditions in the observation hut 
on the other hand could not be established" (see page 175 of ref. [68J). 

Most surprisingly, however, Shankland et al. seem not to realize that Miller's average value 
^^xp ^ 0.044, obtained after their own re-analysis of his observations at Mt. Wilson, when 
compared with the reference classical value ^2^^*^*^ = 0.56 for his apparatus, was giving the same 
observable velocity fobs ~ 8.4 km/s obtained from Miller's re-analysis of the Michelson-Morley 
experiment in Cleveland. Conceivably, their emphasis on the role of temperature effects would 
have been re-considered had they realized the perfect identity of two determinations obtained 
in completely different experimental conditions. In this sense, an interpretation in terms of a 
temperature gradient is only acceptable provided this gradient represents a non-local effect, 
as in our model of the ether drift from a fundamental vacuum energy-momentum flow. 

Another criticism of Miller's work was recently presented by Roberts [80]. This author, 
using the set of data reported in Fig. 8 of ref. [65], raises several objections to the validity 
of Miller's observations. The two main objections concern i) the subtraction of the steady 
thermal drift, which was approximated by Miller as a pure linear effect, and ii) the statistical 
significance of the measurements. Concerning remark i), Roberts reports in his Fig. 3 a broken 
line that reproduces the expected linear trend. He also reports some chosen points (differing 
from the corners of the broken line by 180 degrees) that, due to the 2nd-harmonic nature 
of the ether-drift effect, should lie on the line. However, this expectation ignores that, as 
already pointed out for the Michelson-Morley experiment, real measurements contains large 
first-harmonic effects. These only cancel when taking symmetric combinations of data at 
the various angles 9 and ir + 9. As a matter of fact, the autocorrelative methods and fur- 
ther tests applied by the Shankland team over all of Miller's data confirmed the linear drift 
approximation as remarkably good (see their footnote 21 on page 177 of [68J). 

Concerning remark ii), according to Roberts, the experimental uncertainties are so large 
that the observed 2nd-harmonic effect has no statistical significance. To check this point 
we have re-computed ourselves the fringe shifts for the set of 20 turns of the interferometers 
(reported in Fig.8 of ref. [65]) considered by Roberts, by following the same procedure explained 
in Sect. 3. The resulting symmetric combinations of fringe shifts 

2A 

are reported in our Table 6. 

We have then fitted these data by including both 2nd and 4th harmonic terms. Notice 
that, differently from Roberts' analysis, we do not perform any averaging of data obtained 
from different turns of the interferometer. For our global fit, to estimate the accuracy of the 
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Table 6: The symmetric combination of fringe shifts B{9) = — ^ — various 

values of for the set of 20 turns of the interferometer reported in Fig. 8 of ref. 165^ . For our 
global fit, following ref.]68^. the nominal accuracy of each entry has been fixed to ±0.050. 



Turn 


0° 


22.5° 


45° 


67.5° 


90° 


112.5° 


135° 


157.5° 


1 


+0.091 


+0.159 


+0.028 


+0.047 


-0.034 


-0.116 


-0.147 


-0.028 


2 


-0.025 


+0.063 


+0.050 


+0.088 


-0.075 


-0.038 


+0.000 


-0.063 


3 


+0.022 


+0.103 


+0.084 


+0.016 


-0.053 


-0.072 


-0.091 


-0.009 


4 


+0.034 


-0.009 


-0.053 


-0.047 


-0.041 


+0.016 


+0.022 


+0.078 


5 


+0.169 


+0.081 


+0.044 


-0.044 


-0.081 


-0.169 


-0.056 


+0.056 


6 


-0.025 


+0.025 


+0.025 


+0.025 


+0.025 


-0.025 


-0.025 


-0.025 


7 


+0.081 


+0.094 


+0.056 


+0.069 


-0.119 


-0.106 


-0.094 


+0.019 


8 


+0.066 


+0.072 


-0.022 


-0.066 


-0.059 


-0.003 


+0.003 


+0.009 


9 


+0.041 


+0.084 


+0.078 


+0.022 


-0.134 


-0.141 


+0.003 


+0.047 


10 


+0.016 


+0.072 


+0.078 


-0.016 


-0.009 


-0.003 


-0.047 


-0.091 


11 


+0.009 


+0.053 


+0.097 


-0.009 


-0.116 


-0.072 


+0.022 


+0.016 


12 


+0.022 


+0.016 


+0.059 


+0.003 


-0.053 


-0.009 


-0.016 


-0.022 


13 


+0.000 


+0.063 


+0.025 


+0.038 


+0.050 


-0.038 


-0.075 


-0.063 


14 


-0.034 


+0.047 


+0.078 


+0.009 


-0.009 


-0.028 


-0.047 


-0.016 


15 


+0.113 


+0.125 


+0.138 


+0.000 


-0.088 


-0.125 


-0.113 


-0.050 


16 


+0.025 


+0.050 


+0.025 


+0.050 


-0.025 


-0.050 


-0.025 


-0.050 


17 


+0.000 


-0.012 


-0.025 


+0.063 


+0.000 


-0.012 


-0.025 


+0.013 


18 


+0.044 


+0.050 


+0.019 


-0.019 


-0.056 


-0.044 


-0.031 


+0.031 


19 


+0.053 


+0.059 


+0.016 


-0.028 


-0.022 


-0.066 


-0.009 


-0.003 


20 


+0.059 


+0.041 


+0.122 


+0.003 


-0.066 


-0.084 


-0.053 


-0.022 
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various determinations, we have followed ref.[68] and adopted a nominal uncertainty ±0.050 
for each entry of Table 6. From the fit, where the 4th harmonic is completely consistent with 
the background {A^^'^ = 0.004 ± 0.012), we have obtained a chi-square of 130 for 157 degrees 
of freedom and the following values 



= 24° ± 7° 



(73) 



Here errors correspond to the overall boundary Ax^ = +3.67, as appropriate for a 70% C. 
L. in a 3-parameter fit [81]. Notice that, even though the fitted A2 Eq. (j73p is only 20% larger 
than the nominal accuracy ±0.050 of each entry, the data are distributed in such a way to 
produce a 5a evidence for a non-zero 2nd harmonic. 

As for Illingworth's experiment, we have also analyzed the results obtained from the 
individual turns of the interferometer. To this end, we report in Figs. 7 and 8 the plots of 
the azimuth and of 
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Figure 7: The azimuth ( in degrees ) for the 20 individual turns of the interferometer reported 
in Table 6. The average uncertainty of each determination is about ±20°. The band between 
the two horizontal lines corresponds to the global fit 9q = 24° ±7°. Each individual value could 
also be reversed by 180 degrees. 

To conclude our analysis of Miller's experiment, we want to mention that other objections 
to the overall consistency of his solution for the Earth's cosmic motion [65j were raised by 
von Laue [82] and Thirring [83]. Their argument, which concerns the observed displacement 
of the maximum of the fringe pattern averaged over all sidereal times, was also re-proposed 
by Shankland et al. [68j and amounts to the following. 

^This probability content assumes a Gaussian distribution as for typical statistical errors. 
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Figure 8: The 2nd-harmonic amplitude for the 20 individual turns of the interferometer 
reported in Table 6. The average uncertainty of each determination is about ±0.030. The 
band between the two horizontal lines corresponds to the global fit A2 = 0.061 ±0.012. Within 
their errors, these individual values correspond to an observable velocity in the range 4^14 
km/s. 

By assuming relations (j3ip — ()43p and denoting by (...) the daily average of any given 
quantity, one finds, at any angle 9, the daily averaged fringe shift 

{^^) = 2{C{t))cos2e (74) 
A 

since {S{t)) = with 

A c2 4 

The result can then be cast into the form 1 68 



{C{t)) = - \ ^ -^ -(3cos27-l)cos^(A (75) 



{^M) = v'F{j,cP)cos29 (76) 

Therefore, since the latitude (p is a constant and the angular declination 7 is fixed at any 
specific epoch, the daily averaged fringe shifts should all have a common maximum at the 
value 6 = 0. Only the amplitude can be different at different epochs. Instead, in Miller's 
observations the location of the maximum was differently displaced from the meridian (see 
Figs. 25 of ref.[65] and Fig. 3 of ref.[68j). The presence of such effect has always represented a 
problem for the overall consistency of Miller's solution for the Earth's cosmic motion [65]. 

However, in this derivation, one assumes that any physical signal should only exhibit the 
smooth modulations expected from the Earth's rotation. As anticipated in Sect. 2, and dis- 
cussed in connections with the Michelson-Morley and Illingworth experiments, one might be 
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faced with the more general scenario where the two velocity components Vx{t) and Vy{t) in 
Eq. (j44p are not smooth periodic functions but exhibit stochastic behaviour. In this differ- 
ent perspective, combining observations of different days and different epochs becomes more 
delicate and there might be non-trivial deviations from Eq. ()76p . We shall therefore conclude 
our analysis of Miller's experiments by recalling the remarkable consistency of the velocity 
value V ~ 350 km/s (obtained from the 2nd-harmonic amplitude A2^^ ~ 0.044 computed 
by the Shankland team) with those from the Michelson-Morley, Morley- Miller and Kennedy- 
Illingworth experiments. In this sense, this bulk of Miller's work will remain. 

7. Michelson-Pease-Pearson 

Let us further compare with the experiment performed by Michelson, Pease and Pearson 
|84l I85j . They do not report numbers so that we can only quote from the original article 
|85j which reports the outcome of the measurements performed in the most refined version of 
the experiment: " In the final series of experiments, the apparatus was transferred to a well- 
sheltered basement room of the Mount Wilson Laboratory. The length of the light path was 
increased to eighty- five feet, and the results showed that the precautions taken to eliminate 
temperature and pressure disturbances were effective. The results gave no displacement as 
great as one-fiftieth of that to be expected on the supposition of an effect due to a motion of 
the solar system of three hundred kilometers per second". On the other hand, in ref.[83], after 
similar comments on the length of the apparatus and on the precautions taken to eliminate 
the various disturbances, one finds this other statement "The results gave no displacement as 
great as one-fifteenth of that to be expected on the supposition of an effect due to a motion 
of the solar system of three hundred kilometers per second. These results are differences 
between the displacements observed at maximum and at minimum at sidereal times, the 
directions corresponding to Dr. Stromberg's calculations of the supposed velocity of the solar 
system". In the same paper, the authors report that, according to Stromberg's calculations 
" a displacement of 0.017 of the distance between fringes should have been observed at the 
proper sidereal times" . 

Clearly, although not explicitly stated, they were assuming that some unknown mechanism 
was largely reducing the fringe shifts with respect to the naive non-relativistic value associated 
with a kinematical velocity of 300 km/s. Thus one could try to conclude that their experiment 
implies fringe shifts < 0.017 ~ 0.001. However this is not what they say (they speak 
of differences between fringe displacements) and, in any case, this interpretation does not fit 
with the result reported by Shankland et al. [6^ (see their Table I). According to these other 
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authors, the typical observed fringe shifts observed by Michelson, Pease and Pearson were of 
the order of ±0.005. 

To try to understand this intricate issue, we have been looking at another article [86] 
which, surprisingly, was signed by F. G. Pease alone. Here, one discovers that, in the first 
stage of the experiment, the fringe shifts had a typical magnitude of about ±0.030. Later on, 
however, by reducing substantially the rotation speed of the apparatus, the observed effects 
became considerably smaller. 

Pease declares that, in their experiment, to test Miller's claims, they concentrated on a 
purely 'differential' type of measurement. For this reason, he only reports the difference 

el^J - ( )5.30 - \ )l7.30 [I') 

between the mean fringe shifts (■^)5.305 obtained after averaging over a large set of obser- 
vations performed at sidereal time 5.30, and the mean fringe shifts (^)i7.3o obtained after 
averaging in the same period at sidereal time 17.30. The quantity e{6) has typical magnitude 
of ±0.004 or smaller. However, as already anticipated in Sect. 2, by averaging observations 
performed at a given sidereal time one is assuming the smooth modulations of the signal 
described by Eqs. (j37p . (j38p . Otherwise, one will introduce uncontrolled errors. For instance 
if, consistently with Illingworth's and Miller's data, there were substantial stochastic compo- 
nents in the signal, the cancelations introduced by a naive averaging process would become 
stronger and stronger by increasing the number of observations. 

Therefore, from these values, nothing can be said about the magnitude of the fringe shifts 
^^p^ obtained, before any averaging procedure and before any subtraction, in individual 
measurements at various hours of the day. Pease reports a plot of just a single observation, 
performed when the length of the optical path was still 55 feet, where the even fringe shift 
combinations Eq. (j50p vary approximately in the range ±0.007. This is equivalent to fringe 
shifts of about ±0.011 with a length of 85 feet and could hardly be taken as indicative 
of the whole sample of measurements. In this situation, one can only adopt the estimate 
A2 ~ 0.010 ± ... for the value of the 2nd-harmonic amplitude, for optical path L=85 feet, 
whose uncertainty cannot be estimated in the absence of information on the other individual 
sessions. Then, for this configuration, where j (30km/s) ^ 0.45, this is equivalent to 

Voi,s = (4.5 ± ...) km/s (78) 

or, by using Eq. ([30|) . to 

= (185 ± ...) km/s (79) 
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We emphasize that Miher's extensive observations, as reported in Fig. 22 of ref . |65] (see also 
our Fig. 8), gave fluctuations of the observable velocity lying, within the errors, in the range 
4—14 km/s which has been smoothed in our Fig. 2. For this reason, even though Miller's 
reconstruction of the Earth's cosmic motion is not internally consistent, a single observation 
which gives t'obs ~ 4.5 km/s does not represent a refutation of the whole Miller experiment. 
This becomes even more true by noticing that the single session selected by Pease, within a 
period of several months, was chosen to represent an example of extremely small ether-drift 
effect. 

8. Joos 

One more classical experiment, performed by Georg Joos in 1930, has finally to be considered. 
For the accuracy of the measurements (data collected at steps of 1 hour to cover the full 
sidereal day that were recorded by photocamera), this experiment cannot be compared with 
the other experiments (e.g. Michelson-Morley, Illingworth) where only observations at few 
selected hours were performed and for which, in view of the strong fluctuations of the azimuth, 
one can just quote the average magnitude of the observed velocity. Moreover, differently from 
Miller's, the amplitudes of all basic Joos' observations can be reconstructed from the published 
articles [871 EH]- As such, this experiment deserves a more refined analysis and will play a 
central role in our work. 

Joos' optical system was enclosed in a hermetic housing and, traditionally, it was always 
assumed that the fringe shifts were recorded in a partial vacuum. This is supported by several 
elements. For instance, when discussing his device for electromagnetic fine movements of the 
mirrors, Joos explicitly refers to the condition of an evacuated apparatus, see p. 393 of |87j . 
This aspect is also confirmed by Miller who, quoting Joos' experiment, explicitly refers to an 
"evacuated metal housing" in his article [65] of 1933. This is particularly important since later 
on, in 1934, Miller and Joos had a public letter exchange [89] and Joos did not correct Miller's 
statement. On the other hand, Swenson [90] explicitly reports that fringe shifts were finally 
recorded with optical paths placed in a helium bath. In spite of the fact that this important 
aspect is never mentioned in Joos' papers, we shall follow Swenson and assume that during 
the measurements the interferometer was filled by gaseous helium at atmospheric pressure. 

The observations were performed in Jena in 1930 starting at 2 P.M. of May 10th and 
ending at 1 P.M. of May 11th. Two measurements, the 1st and the 5th, were finally deleted 
by Joos with the motivation that there were spurious disturbances. The data were combined 
symmetrically, in order to eliminate the presence of odd harmonics, and the magnitude of the 
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Figure 9: The selected set of data reported by Joos 187[ \88^ . The yardstick corresponds to 
1/1000 of a wavelength so that the experimental dots have a size of about 0.4 • 10^^. This 
corresponds to an uncertainty ±0.2 • 10^'^ in the extraction of the fringe shifts. 

fringe shifts was typically of the order of a few thousandths of a wavelength. To this end, 
one can look at Fig. 8 of [88] (reported here as our Fig. 9) and compare with the shown size of 
1/1000 of a wavelength. From this picture, Joos decided to adopt 1/1000 of a wavelength as 
an upper limit and deduced an observable velocity Vohs ^ 1-5 km/s. To derive this value, he 
used the fact that, for his apparatus, an observable velocity of 30 km/s would have produced 
a 2nd-harmonic amplitude of 0.375 wavelengths. 

Still, since it is apparent from Fig. 9 that some fringe displacements were definitely larger 
than 1/1000 of a wavelength, we have decided to extract the values of the 2nd-harmonic 
amplitude A2 from the 22 pictures. Differently from the values of the azimuth, this can be 
done unambiguously. The point is that, due to the camera effect, it is not clear how to fix the 
reference angular values in Fig. 9 for the fringe shifts. Thus, one could choose for instance the 
set (k=l, 2, 3, 4) 9k =(0°, 45°, 90°, 135°) or the different set Ok =(360°, 315°, 270°, 225°). 
Or, by noticing that in Fig. 9 there is a small misalignment angle 6* ~ 17° (which actually 
from [87J might instead be 22.5°) between the dots of Joos' fringe shifts and the N, W, and S 
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marks, one could also adopt other two set of values, namely 9k ={0° + 0* , 45° + 9*, 90° + 9*, 
135° + 9*) or 9k ={360° - 9*, 315° - 9*, 270° - 9*, 225° - 9*). By fitting the fringe shifts 
of Fig. 9 to the 2nd-harmonic form Eq. (j56p . these four options for the reference angles 9^, 
would give exactly the same amplitude A2 but four different choices for the azimuth, i.e. —9o, 
—9q + 0*, 9q — 9* and ^o- This basic ambiguity should be added to the standard uncertainty 
in the azimuth that, due to the 2nd-harmonic nature of the measurements, could always be 
changed by adding ±180 degrees Therefore, since clearly there is only one correct choice 

for the angles 9^^ '^'^^^ \\i^\7(^ Tirt^^t^rrt^A nnt tn minfp tVipnrpfir'pl nnr-prfpintipc: on tVi^^ azimuth and 

just concentrate o 1 Fig. 10. 
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Figure 10: Joos' 2nd-harmonic amplitudes, in units 10 The vertical band between the two 
lines corresponds to the range (1.4 it 0.8) • 10~^. 

By computing mean and variance of the individual values, we obtain an average 2nd- 
harmonic amplitude 

(4°°') = (1.4 ± 0.8) • 10"=^ (80) 
and a corresponding observable velocity 

^^obs ~ 1.8t[J| km/s (81) 

By correcting with the helium refractive index, Eqs. ()30p and (I8ip would then imply a true 
kinematical velocity v ~ 217^^9 km/s. 

*As an example, one can consider the azimuth for Joos' picture 20. Depending on the choice of the reference 
angles 9k, one finds 0o ~ 329°, 329° + 9*, 31° - , 31° or 9o ~ 149°, 149° + 9*, 211° - 9", 211°. 
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Tabic 7: The 2nd-harmonic amplitude obtained from the 22 Joos pictures of our Fig. 9. The 
uncertainty in the extraction of these values is about ±0.2-10^^ (the size of the dots in Fig. 9). 
The mean amplitude over the 22 determinations is (^2°°*^) = 1-4 • 10~^. 
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However, this is only a first and very partial view of Joos' experiment. In fact, we have 
compared Joos' amplitudes with theoretical models of cosmic motion. To this end, one has 
first to transform the civil times of Joos' measurements into sidereal times. For the longitude 
11.60 degrees of Jena, one finds that Joos' observations correspond to a complete round in 
sidereal time in which the value t = 0° = 360° is very close to Joos' picture 20. Then, by 
using Eqs. ([3T]) and (jMj) . one can use this input and compare with theoretical predictions for 
the amplitude which, for the given latitude (p = 50.94 degrees of Jena, depend on the right 
ascension a and the angular declination 7. To this end, it is convenient to first re-write the 
theoretical forms as 

A,(t) co= 2S„(t) = 2C(t) = .A>_jA> . 2.6 . 10- ^^^^ (82) 

and 



Mt)s,.mt) = 2S(*) = . 2.6 . 10-3 (83) 



where we have used the numerical relation for Joos's experiment j (30km/b) ^ 0.375 and the 
value of the helium refractive index. Then, by approximating Vx{t) ~ Vx{t), Vy{t) ~ Vy{t) and 



using Eg. ([3^ for the scalar combination v{t) = y 'ul(t) + we have fitted the amplitude 

data of Table 7 to the smooth form 

Ar°°*^(i) = const • sin^ z{t) (84) 

where cosz(t) is defined in Eq. (j3ip . The results of the fit[^ 

a = 168° ±30° 7 = -13° ±14° (85) 

confirm that, as found in connection with the Illingworth experiment, the Earth's motion 
with respect to the CMB (which has a ~ 168° and 7 ~ —6°) could serve as a useful model to 
describe the ether-drift data. 

Still, in spite of the good agreement with the CMB a— and 7— values obtained from the 
fit Eq. ()85p . the nature of the strong fluctuations in Fig. 10 remains unclear. Apart from this, 
there is also a sizeable discrepancy in the absolute normalization of the amplitude. In fact, 
by assuming the standard picture of smooth time modulations, the mean amplitude over all 
sidereal times can trivially be obtained from the mean squared velocity Eq.(| 



^Actually, there is another degenerate minimum at a = 348° ± 30° and 7 = 13° ± 14° because sin^ z(t) 
remains invariant under the simultaneous replacements a a-\- 180° and 7 — ^ —7. However, due to the close 
agreement with the CMB parameters we have concentrated on solution H85p . 
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For the CMB and Jena, this gives \J [v^) ~ 330 km/s so that one would naively predict from 
Eqs.dSS]), dSSD 

(A— (t)) ~ 2.6 • 10-3 ^3^^^^^^]^^, ~ 3.2 . 10-3 (87) 

to be compared with Joos' mean value (^2°°'') — (l-4ib0.8)-10-3. In the standard picture, this 
experimental value leads to the previous estimate \J [v^) 217 km/s and not to ~ 330 

km/s so that it is necessary to change the theoretical model to try to make Joos' experiment 
completely consistent with the Earth's motion with respect to the CMB. 

To try to solve this problem, and understand the origin of the observed strong fluctuations, 
we have used the same model Eqs. (f60|) . ([6T]) of Sect. 5, to simulate stochastic variations of the 
velocity field. As anticipated however, in this case of the Joos experiment, we have modified 
the theoretical framework. Namely, we have allowed the two random parameters x„(l) and 
Xn(2) to vary in the range [—Vx{t),Vx(t)] and the other two parameters yn(l) and yn(2) to 
vary in the different range [—Vy{t),Vy{t)], where Vx{t) and Vy{t) are defined in Eqs. ([3T]) — (p3|) . 
In this way, for each time t, Eqs. (|62p now become 

{xl{t = 1, 2))stat = {yl[l = 1, 2))stat = ^ (88) 

It is understood that the latitude corresponds to Joos' experiment while V , a and 7 describe 
the Earth's motion with respect to the CMB. Notice that, in this model, there will be a 
substantial reduction of the amplitude with respect to its smooth prediction. To estimate 
the order of magnitude of the reduction, one can perform a full statistical average (as for an 
infinite number of measurements) and use Eqs. (j88p in Eqs. (j82p . (j83p for our case rj = \. This 
gives 

(^2(t))stat ~ 2.6 • 10-3 - — —^-^ - S^ — = — Af^^^'^it) (89) 
^ (300 km/s) 2 3^n2 18 ^ ^' ^ ' 

By also averaging over all sidereal times, for the CMB and Jena, one would now predict a 
mean amplitude of about 1.7 • 10^^ and not of 3.2 • lO^^. 

After having fixed all theoretical inputs, we have analyzed the dependence of the numerical 
results on the remaining parameters of the simulation, namely the number of Fourier modes 
(in the available range A^ < 10^) and the integer number s (the 'seed') which determines the 
random sequence. In particular, the dependence on the latter is usually quoted as theoretical 
uncertainty. For this reason, for Illingworth's experiment in Sect. 5 we had produced several 
copies of the high-statistics simulation in Fig. 6 by quoting values for the standard deviations 
Eq. (j6l|) which take into account the observed s— dependence of the results. 
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Table 8: The 2nd-harmonic amplitude obtained from a single simulation of 22 instantaneous 
measurements performed at Joos ' times. The stochastic velocity components are controlled by 
the kinematical parameters {V, a, 7)cmb o,s explained in the text. The mean amplitude over 
the 22 determinations is (^f™"') = 1.38 • lO^^. 
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Here, we have started by doing something similar. However, since it is not possible to 
consider at once all characteristics of a given configuration, we have first concentrated on 
the simplest statistical indicator, namely the mean amplitude (Ag™"^) obtained by averaging 
over all sidereal times. Quite in general, this can be evaluated for a variety of configurations 
which depend on the number n of measurements that one wants to simulate and the interval 
At between two consecutive measurements. For instance, Joos' experiment corresponds to 
n = 24 (actually n = 22 since Joos finally deleted two observations) and At ~ 3600 seconds. 
At the same time, the simulations become quite lengthy for large N, large n and small At. 
Therefore, we have first performed a scan of s— values for N = 10^ and then studied a few s 
by increasing N. To give an idea of the spread of the central values, due to changes of the 
pair {N,s), we report below the approximate results of this analysis for some choices of the 
pair (n. At) 

(^|™^^(n = 24, At = 3600 s)) ~ (1.7 ± 0.8) • lO'^ (90) 

(^|™"'(n = 1440, At = 60 s)) ~ (1.7 ± 0.3) • lO'^ (91) 

(yl|'"^^'(n = 240, At = 3600 s)) ~ (1.8 ± 0.5) • 10"^ (92) 

As it might be expected, the average (Ag™"^) becomes more stable by increasing the number 
of observations. Concerning the individual values ^2™"'(^j)) with i = l,..,n, they have a 
large spread, about (1-^4) • 10^'^. This is in agreement with the 'fat-tailed' distributions of 
instantaneous values expected in turbulent flows [TBI [77| (compare with Fig. 6 in Sect. 5). 
However this other spread can be reduced by starting to average the data in some interval 
of time to. In this case, the spread of the resulting average values (^2 to decreases as 
We emphasize that, by performing extensive simulations, there are occasionally very 
large spikes of the amplitude at some sidereal times, of the order (10 20) • 10^'^. The effect 
of these spikes gets smoothed when averaging over many configurations but their presence is 
characteristic of a stochastic-ether model. With a standard attitude, where the ether drift 
is only expected to exhibit smooth time modulations, the observation of such effects would 
naturally be interpreted as a spurious disturbance (Joos' omitted observations 1 and 5?). 

After this preliminary study, we have then concentrated on the real goal of our simulation, 
i.e. to compare with the single Joos configuration of 22 entries in Table 7. To this end, one 
could first try to look for the 'best seed', or subset of seeds, which can minimize the difference 
between the generated configurations and Joos' data. This standard task, usually accom- 
plished by minimizing a chi-square, is difficult to implement here. In fact, it is problematic to 
construct a function x^(s) and look for its minima because a seed s and the closest seeds s±l 
give often vastly different configurations and chi-square. For this reason, we have followed 
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Figure 11: Joos' experimental amplitudes in Table 7 are compared with the single simulation 
of 22 measurements for fixed {N,s) in Table 8. By changing the pair (N,s), the typical 
variation of each simulated entry is (1-^4)- 10^^ depending on the sidereal time. We also 
show two 5th-order polynomial fits to the two different sets of values. 

an empirical procedure by forming a grid and selecting a set of seeds whose mean amplitude 
(for n = 24 and At = 3600 s) gets close to Joos's mean amplitude (^2°°**) = 1-4 • 10~^ for a 
large number N of Fourier modes. One of such seeds gave a sequence (Af™"') =1-66, 1.40, 

I. 08, 1.21 and 1.38 (in units lO'^), for N = 10^ 10^, 10^ 10^ and 5 • 10^ respectively, and 
the configuration with = 5 • 10^ was finally chosen to give an idea of the agreement one 
can achieve between data and a single numerical simulation for fixed {N,s). The simulated 
values are reported in Table 8 and a graphical comparison with Joos' data is shown in Fig. 

II. We emphasize that one should not compare each individual entry with the corresponding 
data since, by changing {N,s), the simulated instantaneous values vary typically of about 
(1-^4) • 10~^ depending on the sidereal time. Instead, one should compare the overall trend of 
data and simulation. To this end, we show two 5th-order polynomial fits to the two different 
sets of values. 

A more conventional comparison with the data consists in quoting for the various 22 
entries simulated average values and uncertainties. To this end, we have considered the mean 
amplitudes (^|™^'(ti)) defined by averaging, for each Joos' time tj, over 10 hypothetical 
measurements performed on 10 consecutive days. For each ti, the observed effect of varying 
(N, s) has been summarized into a central value and a symmetric error. The values are 
reported in Table 9 and the comparison with Joos' amplitudes is shown in Fig. 12. 
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Table 9: The 2nd-harmonic amplitudes obtained by simulating the averaging process over 
10 hypothetical measurements performed, at each Joos' time, on 10 consecutive days. The 
stochastic velocity components are controlled by the kinematical parameters {V, a, 7)cmb as 
explained in the text. The effect of varying the pair (iV, s) has been approximated into a 
central value and a symmetric error. The mean amplitude over the 22 determinations is 
^^simui^ = 1.8- 10-^ 
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Figure 12: Joos' experimental amplitudes in Table 7 are compared with our simulation in 
Table 9. 

The spread of the various entries is larger at the sidereal times where the projection at Jena 
of the cosmic Earth's velocity becomes larger. The tendency of Joos' data to lie in the lower 
part of the theoretical predictions in Table 9 mostly depends on our use of symmetric errors. 
In fact, by comparing in some case with the histograms of the basic generated configurations 
^2™"'(^j)i have seen that our sampling method of (^2™^'(^i))> based on a grid of {N,s) 
values, typically underestimates the weight of the low-amplitude region in a prediction at the 
70% C.L. . This can also be checked by considering the single simulation of Table 8 and 
counting the sizeable fraction of amplitudes ^2™^'(tj) < 0.5 • 10~^. For this reason, one could 
improve the evaluation of the probability content. However, in view of the good agreement 
already found in Fig. 12 (x^ = 13/22), we did not attempt to carry out this more refined 
analysis. 

In conclusion, after the first indication obtained from the fit Eq. (j85|) . we believe that the 
link between Joos' data and the Earth's motion with respect to the CMB gets reinforced by 
our simulations. In fact, by inspection of Figs. 11 and 12, the values of the amplitudes and the 
characteristic scatter of the data are correctly reproduced. In principle, there could be space 
for further refinements by taking into account the Earth's orbital motion in the input values 
for V, a and 7. 

From this agreement with the data, we then deduce that the previous value for the kine- 
matical velocity v ~ 217]i^7g km/s, obtained by simply correcting with the helium refractive 
index the average observable velocity ()8ip . has to be considerably increased if one allows for 
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stochastic variations of the velocity field. In fact, the magnitude of the fluctuations in Vx and 



Vy is controlled by the same scalar parameter v{t) = yv'^{t) + v'^{t) of Eq.([33D- In view of 
the good agreement between data and our numerical simulations, we conclude that Joos' data 
are consistent with a range of kinematical velocity v = SSO^^q km/s which corresponds to 
Eq.dMD for 4> = 50.94°, V = 370 km/s, a = 168° and 7 = -6°. 

9. Summary and conclusions 

The condensation of elementary quanta and their macroscopic occupation of the same quan- 
tum state is the essential ingredient of the degenerate vacuum of present-day elementary 
particle physics. In this description, one introduces implicitly a reference frame S, where the 
condensing quanta have k = 0, which characterizes the physically realized form of relativ- 
ity and could play the role of preferred reference frame in a modern re-formulation of the 
Lorentzian approach. To this end, we have given in the Introduction some general theoretical 
arguments related to the problematic notions of a non-zero vacuum energy and of an exact 
Lorentz-invariant vacuum state. These arguments suggest the possibility of a tiny vacuum 
energy-momentum flux, associated with an Earth's absolute velocity v, which could affect 
differently the various forms of matter. Namely, it could produce small convective currents 
in a loosely bound system such as a gas or dissipate mainly by heat conduction with no ap- 
preciable particle flow in strongly bound systems as liquid or solid transparent media. In the 
former case, by introducing the refractive index M of the gas, convective currents of the gas 
molecules would produce a small anisotropy, proportional to {M — l)(v/c)^, of the two-way 
velocity of light in agreement with the general structure Eq. (j23p or with its particular limit 
Eq. (|24p . Notice that this tiny anisotropy refers to the system S' where the container of the 
gas is at rest. In this sense, contrary to standard Special Relativity, S' might not define a 
true frame of rest. This conceptual possibility can be objectively tested with a new series 
of dedicated ether-drift experiments where two orthogonal optical resonators are filled with 
various gaseous media by measuring the fractional frequency shift ^vjv between the two 
resonators. By assuming the typical value v ~ 300 km/s of most cosmic motions, one expects 
frequency shifts < 10^^'^ for gaseous helium and < 10^^ for air, which are well within the 
present technology. 

Given the heuristic nature of our approach, and to further motivate the new series of 
dedicated experiments, we have tried to get a first consistency check. In fact, by adopting 
Eq. (|24p . the frequency shift between the optical resonators is governed by the same classical 
formula for the fringe shifts in the old ether-drift experiments with the only replacement 
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— )• 2{M —X)v^ . Thus testing our picture is very simple: one should just check the consistency 
of the kinematical f's obtained in different experiments. 

In the old times, experiments were performed with interferometers where light was propa- 
gating in gaseous media, air or helium at atmospheric pressure, where [M — 1) is a very small 
number. In this regime, the theoretical fringe shifts expected on the basis of Eqs. ()23p and 
(|24p are much smaller than the classical prediction {v/cf' . Another important aspect of these 
classical experiments is that one was always expecting smooth sinusoidal modulations of the 
data due to the Earth's rotation, see Eqs. ([35]) . ([37|) and ([38]) . As emphasized in Sect. 2, we 
now understand the logical gap missed so far. The relation between the macroscopic Earth's 
motions (daily rotation, annual orbital revolution,...) and the ether-drift experiments depends 
on the physical nature of the vacuum. Assuming Eqs. (|37p and (j38p . to describe the effect of 
the Earth's daily rotation, amounts to considering the vacuum as some kind of fluid in a 
state of regular, laminar motion for which global and local properties of the flow coincide. 
Instead, several theoretical arguments (see e.g. refs.p ^ [T5 t [T6 l [T7 1 [5611571 [58 } [59 1 160]) suggest 
that the physical vacuum might behave as a stochastic medium similar to a turbulent fluid 
where large-scale and small-scale motions are only indirectly related. In this case, there might 
be non-trivial implications. For instance, due to the irregular behaviour of turbulent flows, 
vectorial observables collected at the same sidereal time might average to zero. However, this 
does not mean that there is no ether-drift. More generally, the relevant Earth's motion with 
respect to S might well correspond to that indicated by the anisotropy of the CMB, but it be- 
comes non trivial to reconstruct the kinematical parameters from microscopic measurements 
of the velocity of light in a laboratory. These arguments make more and more plausible that 
a genuine physical phenomenon, much smaller than expected and characterized by stochastic 
variations, might have been erroneously interpreted as an instrumental artifact thus leading 
to the standard 'null interpretation' of the experiments reported in all textbooks. 

Now, our analysis of Sects. 3— 8 shows that this traditional interpretation is far from obvi- 
ous. In fact, by using Eqs. fp^ . (j29l) and ([30]) . the small residuals point to an average velocity 
of about 300 km/s, as in most cosmic motions. In this alternative interpretation, the indica- 
tions of the various experiments are summarized in our Table 10 0. As a summary of our 

Other determinations of less accuracy could also be included, as for the 1881 Michelson experiment in 
Potsdam [91] or Tomaschek's starlight experiment [92] or the Piccard and Stahel experiment which was first 
performed in a ballon [93] and later [94] on the summit of Mt. Rigi in Switzerland. These results were 
summarized in Table I of ref . [68j and by Miller [65]. In the 1881 Potsdam experiment the fringe shifts were 
in the range 0.002 0.007 to be compared with an expected 2nd-harmonic of 0.02 for 30 km/s. This means 
observable velocities (9 ~ 18) km/s which are comparable and even larger than those of the 1887 experiment. 
In Tomaschek's starlight experiment, fringe shifts were about 15 times smaller than those classically expected 
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Table 10: The average velocity observed (or the limits placed) by the classical ether-drift 
experiments in the alternative interpretation of Eqs. (f^^, i2y\} . 



Experiment 


gas in the interferometer 


fobs(km/s) 


t;(km/s) 


Michelson-Morley(1887) 


air 


8.4l};f 




Morley-Miller(1902-1905) 


air 


8.5 ± 1.5 


353 ± 62 


Kennedy(1926) 


helium 


< 5 


< 600 


Illingworth(1927) 


helium 


3.1 ± 1.0 


370 ± 120 


Miller(1925-1926) 


air 




040+79 


Michelson- Pease- Pearson( 1929) 


air 


4.5± ... 


185 ± ... 


Joos(1930) 


helium 


1 8+0-5 


3301^° 



work, we emphasize the following points: 

i) an analysis of the individual sessions of the original Michelson- Morley experiment, in 
agreement with Hicks |66j and Miller |65j (see our Figs. 1 and 2), gives no justification to its 
standard null interpretation. As discussed in Sect. 3, this type of analysis is more reliable. In 
fact, averaging directly the fringe displacements of different sessions requires two additional 
assumptions, on the nature of the ether-drift as a smooth periodic effect and on the absence 
of systematic errors introduced by the re-adjustment of the mirrors on consecutive days, that 
in the end may turn out to be wrong. 

ii) one gets consistent indications from the Michelson-Morley, Morley-Miller, Miller and 
Illingworth-Kennedy experiments. In view of this consistency, an interpretation of Miller's 
observations in terms of a temperature gradient |68j is only acceptable provided this gradient 
represents a non-local effect as in our picture where the ether-drift is the consequence of a 
fundamental vacuum energy-momentum flow. We have also produced numerical simulations 
of the Illingworth experiment in a simple statistically isotropic and homogeneous turbulent- 
ether model. This represents a zeroth-order approximation and is useful to illustrate basic 
phenomenological features associated with the picture of the vacuum as an underlying stochas- 
tic medium. In this scheme, Illingworth's data are consistent with fluctuations of the velocity 
field whose absolute scale is controlled hy v = Vcmb ~370 km/s, the velocity of the Earth's 

for an Earth's velocity of 30 km/s. This gives Vohs ^ 7.7 km/s or v < 320 km/s. From Piccard and Stahel, 
in the most refined version of Mt. Rigi, one gets an observable velocity Vohs ^ 1-5 km/s. Since their optical 
paths were enclosed in an evacuated enclosure, this very low value can easily be reconciled with the typical 
kinematical velocity v ~ 300 km/s of the most accurate experiments in Table 10. 
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motion with respect to the CMB. 

iii) on the other hand, there is some discrepancy with the experiment performed by 
Michelson, Pease and Pearson (MPP). However, as discussed in Sect. 7, the uncertainty cannot 
be easily estimated since only a single basic MPP observation is explicitly reported in the 
literature. Therefore, since Miller's extensive observations (see Fig. 22 of ref.|65j and our 
Fig. 8), within their errors, gave fluctuations of the observable velocity in the wide range 4—14 
km/s, a single observation giving Uobs ~ 4.5 km/s cannot be interpreted as a refutation. This 
becomes even more true by noticing that the single session selected by Pease, within a period 
of several months, was chosen to represent an example of extremely small ether-drift effect. 

iv) some more details are needed to account for the Joos observations. This experiment is 
particularly important since the data were collected at steps of 1 hour to cover the full sidereal 
day and were recorded by photocamera. For this reason, Joos' experiment is not comparable 
with other experiments (e.g. Michelson-Morley, Illingworth) where only observations at few 
selected hours were performed and for which, in view of the strong fluctuations of the azimuth, 
one can just quote the average magnitude of the observed velocity. Moreover, differently from 
Miller's, the amplitudes of all Joos's observations can be reconstructed from the published 
articles \87\ I88j . For these reasons, this experiment has deserved a more refined analysis 
and is central for our work. As discussed in Sect. 8, due to uncertainties in the original 
data analysis, the standard 1.5 km/s velocity value quoted for this experiment should be 
understood as an order of magnitude estimate and not as a true upper limit. Instead, our 
reported observable velocity Vohs ~ l-^log km/s has been obtained from a direct analysis 
of Joos' fringe shifts. From this value, to deduce a kinematical velocity, one still needs the 
refractive index. The traditional view, motivated by Miller's review article [65J and Joos's 
own statements in ref . |87j . is that the experiment was performed in an evacuated housing. In 
these conditions, it would be easy to reconcile a large kinematical velocity v ~ 350 km/s with 
the very small values of the observable velocity. On the other hand, Swenson [90] explicitly 
reports that fringe shifts were finally recorded with optical paths placed in a helium bath. 
Since Joos' papers do not provide any definite clue on this aspect, we have decided to follow 
Swenson's indications. In this case, by simply correcting with the helium refractive index the 
result fobs ~ l-Sli^o'g km/s, one would get a kinematical velocity v ~ 217^yg km/s. However, 
as discussed in detail in Sect. 8, this is only a first partial view of Joos' experiment. In fact, 
by fitting the experimental amplitudes in Table 7 to various forms of cosmic motion (see 
Eq. (j85p ) we have obtained angular parameters which are very close to those that describe 
the CMB anisotropy (right ascension acMB ~ 168° and angular declination 7cmb ~ —6°). 
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Still, to get a complete agreement, one should explain the absolute normalization of the 
amplitudes and the strong fluctuations of the data. Thus we have improved our analysis 
by performing various numerical simulations where the velocity components in the plane of 
the interferometer Vx{t) and Vy{t), which determine the basic functions C{t) and S{t) through 
Eqs. ()44p and the fringe shifts through Eq. (j35p . are not smooth functions but are represented as 
turbulent fluctuations. Their Fourier components in Eqs. (j60p and (j6ip now vary within time- 
dependent ranges Eqs. ([52|) — ([55]) . [—Vx{t),Vx{t)] and [—Vy{t),Vy{t)] respectively, controlled by 
the macroscopic parameters (V^, a,7)cMB- Taking into account these stochastic fluctuations 
of the velocity field tends to increase the fitted average Earth's velocity, see Eq. (j89p . and can 
reproduce correctly Joos' 2nd-harmonic amplitudes and the characteristic scatter of the data, 
see Figs. 11 and 12. In view of this consistency, we conclude that the range v = SSOl^g km/s 
(corresponding to Eq. (j34p for CMB and Joos' laboratory) is actually the most appropriate 
one. 

The more refined analysis adopted for the Joos experiment provides an explicit example 
of the previously mentioned non-trivial ingredients that might be required to reconstruct the 
global Earth's motion from microscopic measurements performed in a laboratory. For this 
reason, the results reported in Table 10, besides providing an impressive evidence for a light 
anisotropy proportional to {M — l)(w/c)^, with the realistic velocity values v ~ 300 km/s of 
most cosmic motions, could also represent the first experimental indication for the Earth's 
motion with respect to the CMB. Due to the importance of this result, and to provide the 
reader with all elements of the analysis, we present in a second Appendix a brief numerical 
simulation of one noon session of the Michelson-Morley experiment. This has been performed 
in the same framework adopted for the Joos experiment where the velocity components Vx{t) 
and Vy{t) in the plane of the interferometer are represented as turbulent fluctuations varying 
within time-dependent ranges controlled by the macroscopic parameters (V^, a,7)cMB- We 
postpone to a future publication the non-trivial task of performing a complete numerical 
simulation of the whole Michelson-Morley experiment and of the Illingworth experiment (with 
its 32 sessions and the associated sets of 20 rotations for each session) where we'll also compare 
the various theoretical schemes mentioned in Sect. 5 to handle the stochastic components of 
the velocity field. 

We emphasize that the simulation reported in our second Appendix corresponds to a 
single configuration whereas taking into account more and more configurations is essential to 
properly estimate theoretical uncertainties (as for the Joos experiment with the results in our 
Table 9 and Fig. 12). Nevertheless, even this very small sample can provide interesting clues 
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on the real data. For instance, the strong scatter of the fringe shifts at the same 0— values in 
consecutive rotations and the good agreement with the experimental azimuths obtained by 
accepting Hicks' interpretation of the observations of July 8th (see Sect. 3). In this sense, this 
brief numerical analysis reinforces the picture of the classical experiments emerging from our 
Table 10. According to the usual view, the theoretical predictions, with a very low velocity 
f ~ 30 km/s, were much larger than the observed values. Instead, in a modern view of the 
vacuum as a stochastic medium, theoretical predictions, for the realistic velocities v ~ 300 
km/s of most cosmic motions, are now well compatible and, sometimes, even smaller than 
the actual outcome of the observations. This latter case simply means that the experimental 
data were also affected by spurious effects such as deformations induced by the rotation of the 
apparatus or local thermal conditions. This gives a strong motivation to repeat these crucial 
measurements with today's much greater accuracy. 

To this end, let us now briefly consider the modern ether-drift experiments. As anticipated, 
in the modern experiments, the test of the isotropy of the velocity of light consists in measuring 
the relative frequency shift Az/ of two orthogonal optical resonators \6b\ I36j. Here, the analog 
of Eq. (|29p . for a hypothetical physical part of the frequency shift (after subtraction of all 
spurious effects), is 



where is the direction of the ether-drift. This can be interpreted within Eq. ijlOSp where 



A/'medium being the refractive index of the gaseous medium filling the optical resonators. Test- 
ing this prediction, requires replacing the high vacuum usually adopted within the optical 
resonators with a gaseous medium and studying the substantially larger frequency shift in- 
troduced with respect to the vacuum experiments. 

As a rough check, a comparison was made [71[5T] with the results obtained by Jaseja et. al 
|95j in 1963 when looking at the frequency shift of two orthogonal He-Ne masers placed on a 
rotating platform. To this end, one has to preliminarily subtract a large systematic effect that 
was present in the data and interpreted by the authors as probably due to magnetostriction 
in the Invar spacers induced by the Earth's magnetic field. As suggested by the same authors, 
this spurious effect, which was only affecting the normalization of the experimental Ai/, can be 
subtracted by looking at the variations of the data. As discussed in refs.[71 El], the measured 
variations of a few kHz are roughly consistent with the refractive index A/hc-No ~ 1.00004 
and the typical variations of an Earth's velocity as in Eq. ()52p . 




(93) 



'medium I ~ ■'^medium 



(94) 
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More recent experiments |44j — [50] have always been performed in a very high vacuum 
where, as emphasized in the Introduction, the differences between Special Relativity and 
the Lorentzian interpretation are at the limit of visibility. In fact, in a perfect vacuum by 
definition A/Vacuum = 1 so that iSvacuum will vanish. Thus one should switch to the new 
generation of dedicated ether-drift experiments in gaseous systems. Our conclusion is that 
these new experiments should just confirm Joos' remarkable observations of eighty years ago. 
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Appendix I 



To derive Eq. ()24p . one should start from Eq. ()23p which describes hght propagation in a 
gaseous system in the presence of convective currents of the gas molecules originating from 
a fundamental vacuum energy-momentum flow. Due to these convective currents, isotropy 
of the velocity of light would only hold if the solid container of the gas and the observer 
were both at rest in the particular reference frame T, where the macroscopic condensation of 
quanta correspond to the state k = 0. This introduces obvious differences with respect to 
the standard analysis. For instance, let us compare with Jauch and Watson [96j who worked 
out the quantization of the electromagnetic field in a moving medium of refractive index J\f. 
They noticed that the procedure introduces unavoidably a preferred frame, the one where the 
photon energy does not depend on the direction of propagation. Their conclusion, that this 
frame is "usually taken as the system for which the medium is at rest", reflects however the 
point of view of Special Relativity with no preferred frame. Instead, one could consider a 
different scenario where, at least in some limit, the angle-independence of the photon energy 
might only hold for some special frame S. 

To discuss this different case, let us first consider a dielectric medium of refractive index 
M whose container is at rest in S. For an observer at rest in this reference frame, light 
propagation within the medium is isotropic and described by 

TT^TTu-f'"' = (95) 

where 

7'^'' = diag(AA2, -1,-1,-1) (96) 

and vr^ denotes the light 4-momentum vector for the S observer. Let us now consider that 
the container of the medium is moving with some velocity V with respect to S and is at rest 
in some other frame S'. By analogy, light propagation within the medium for the observer in 
S' will be described by 

P^.Pu9^'' = 0, (97) 

where = {E/c,p) and g^'^ denote respectively the light 4-momentum and the effective 
metric for S' . On this basis, by introducing the 5' dimensionless velocity 4- vector = 
{vP,\ /c) (with u^u^ = 1), one can define a transformation matrix A^y = A^y{u^,M) and 
express 

g^.u ^ A^„A\-i-P (98) 
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In this context, requiring the consistency of vacuum condensation with Special Relativity 
corresponds to place all reference frames on the same footing and assume g'^'^ = or 



(99) 



This identification is independent of the physical nature of the medium, being valid for gaseous 
systems as well as for liquid or solid transparent media. In this sense Special Relativity, by 
construction, cannot describe light propagation in the presence of a vacuum energy- momentum 
flow which could affect the various forms of matter differently. 

Instead, consistently with the basic ambiguity in the interpretation of relativity mentioned 
in the Introduction, and with Lorentz' point of view ^ ( "it seems natural not to assume at 
starting that it can never make any difference whether a body moves through the ether or 
not" ) , one could adopt different choices without pretending to determine a priori the outcome 
of any ether-drift experiment. Thus, by noticing that we have at our disposal two matrices, 
namely 5u and the Lorentz transformation matrix A'^i^ associated with V, one could still 
maintain Eq. (j99p for strongly bound systems, such as solid or liquid transparent media, where 
the small energy flux generated by the motion with respect to S should mainly dissipate by 
heat conduction with no appreciable particle flow and no light anisotropy in the rest frame of 
the medium. One could instead identify 



to solve non-trivially the equation g^^'^ = ^^'^ when A/" = 1, i.e. when light propagates 
in vacuum and "f^^ reduces to the Minkowski tensor r]^^'^ . But then, by continuity, it is 
conceivable that Eq. (|lUU|) . up to higher-order terms, can also describe the case N = 1 + e 
of gaseous media. This choice provides a simple interpretative model and a particular form 
of the more general structure Eq. (j23p which corresponds to the Robertson-Mansouri-Sexl 
(RMS) \42\ [33] parametrization for the two-way velocity of light. Moreover, when comparing 
with experiments with optical resonators, the resulting frequency shift is governed by the 
same classical formula for the fringe shifts in the old ether-drift experiments with the only 
replacement V"^ ^ 2{N - l)V'^ . To see this, let us compute g'^^ through the relation 



A'',{u^,N = l) = K 



(100) 



9' 



(101) 



with as in Eq. (j96p . This gives the effective metric for S' 



9' 



(102) 
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with 



1 



In this way, Eq. ([97|) gives a photon energy (uq = 1 + V^/c^ 

^(|P|,^) = 

with 



-KUoC + V|p|2(l + 



1 + KUq 



C = P 



V 

c 



|p|/3 cos ( 



(103) 



(104) 



(105) 



where (3 = '-^ and 9 = ^lab indicates the angle defined, in the laboratory 5' frame, between 
the photon momentum and V. By using the above relation, one gets the one-way velocity of 
hght 

cos 6 + ^/T+~K+~KP^'^I^ 



E( 



IPI 



or to 0{k) and 0{f3^ 



1 + k(1 + /32 



l-K/3cose--/32(l + cos2 9) 



(106) 



(107) 



From this one can compute the two-way velocity 

2c^{9)c^{'K + 1 



-c,{9) 



C-y{9) + €^{71 + 9) 

c 



1-/3^ 



■ sm 



which, as anticipated, is a special form of the more general Eg. ([23 



RMS anisotropy parameter B VA 



{7r/2 + 9) -c^{9) 



~ i3^cos2( 



70 J 



(108) 

We can then define the 
(109) 



There is a subtle difFerence between our Eqs. (|107|) and (|108 |l and the corresponding Eqs. (6) and (10) 
of ref. '7' that has to do with the relativistic aberration of the angles. Namely, in ref. 71, with the (wrong) 
motivation that the anisotropy is 0{l3^), no attention was paid to the precise definition of the angle between 
the Earth's velocity and the direction of the photon momentum. Thus the two-way velocity of hght in the 
S' frame was parameterized in terms of the angle 6 = 6^ as seen in the E frame. This can be explicitly 
checked by replacing in our Eqs. (|107p and p08 |) the aberration relation cos^iab ~ (— /? + cos — /3cos6e) 

or equivalently by replacing cos^e — {13 + cosSiab)/(l + /3cosSiab) in Eqs. (6) and (10) of ref. 7 . However, 
the apparatus is at rest in the laboratory frame, so that the correct orthogonality condition of two optical 
cavities at angles 9 and 71/2 + 9 is expressed in terms of 6 = ^lab and not in terms of 6 — Ot.- This trivial 
remark produces however a non-trivial difference in the value of the anisotropy parameter. In fact, the correct 
resulting \B\ Eq. piO|) is now smaller by a factor of 3 than the one computed in ref.[7] by adopting the wrong 
definition of orthogonality in terms of 6 = 6s- 
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where the pair (v, 9o) describes the projection of V onto the relevant plane and 

1^1 ~ I (AA - 1) (110) 

Prom the previous analysis, by replacing the two-way velocity in Eq. ()28p . one finally obtains 
the observable velocity 

vl^^ ~ 2\B\v'^ ~ 2(Ar - 1)?;^ (Ill) 

to be used in Eq. (|29p . In this way, where one re-obtains the classical formulas with the only 
replacement v — )• Vobs, testing the present scheme requires to check the consistency of the 
kinematical v's obtained in different experiments. 
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Appendix II 



In this second Appendix we'll report the results of a single simulation of an individual 
noon session of the Michelson-Morley experiment. This will be performed within the same 
stochastic-ether model described in Sect. 8 for the Joos experiment. For sake of clarity, we 
recapitulate the various steps so that an interested reader can also run his own simulations. 

One should first express the functions C{t) and S(t) as in Eqs.dM]) and model the two 
velocity components Vx(t) and Vy{t) as in Eqs. ([60|) and (j6T]) . A basic input value is the 
sidereal time of the observation. This has to be inserted, together with the CMB kinematical 
parameters Vcmb ~ 370 km/s, acMB ~ 168°, 7cmb ~ —6°, in Eqs. ([3T]) — (j33|) to fix the 
boundaries [—Vx{t),Vx{t)] and [—Vy{t),Vy{t)] respectively for the random parameters = 
1,2) and = 1,2) entering Eqs. ([5U|) and ([UT]) . In the end, with the simulated C{t) and 
S{t), one should form the fringe shift combination 

= 2C{t) cos 26 + 2S{t) sin 29 = A2{t) cos 2(0 - Ooit)) (112) 

A 

As recalled in Sect. 3, an individual session of the Michelson-Morley experiment consisted of 
6 rotations. Each complete rotation of the interferometer took 6 minutes and the consecutive 
readings of the fringe shifts were performed every 22.5 degrees. Therefore, two consecutive 
readings differed by 22.5 seconds. In these conditions, a numerical simulation of a single 
rotation consists in generating 16 pairs [C{t),S{t)] at steps of 22.5 seconds. 

As the central time of the observations, we have chosen 12 A. M. of July 10, 1887 which, for 
Cleveland, corresponds to a sidereal time r ~ 102°. To select the parameters of the simulation, 
we have compared with the traditional analysis of the experiment where one performs a fit 
to the fringe shifts obtained by averaging the results of the various experimental sessions. In 
our case, averaging the data of the three noon sessions in our Table 1 gives a 2nd-harmonic 
amplitude 

^f*(average data - noon) ~ 0.012 (113) 
We have thus considered the exact amplitude 

^exact(^) = 2^/C^t)+S^t) (114) 

and selected a particular configuration whose global average over the 6 turns gives {A2^^^^{t)) ~ 
0.012. Of course, this condition can be realized by a very large number of configurations. 
These can produce very different fringe shifts at the same 0— values and sizeable variations of 
the fitted amplitude and azimuth. Taking into account these variations is essential to perform 
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a full numerical simulation and estimate theoretical uncertainties (as done for the Joos ex- 
periment with Table 9 and Fig. 12). However, our intention here is just to give an idea of the 
agreement one can achieve between data and a single numerical simulation. We thus postpone 
to a future publication a complete analysis of the whole Michelson-Morley experiment and 
of the Illingworth experiment (with its 32 sessions and the associated sets of 20 rotations for 
each session) where we'll also compare the various theoretical schemes mentioned in Sect. 5 to 
handle the stochastic components of the velocity field. 

The results of our single simulation for [2C{t),2S{t)] are reported in Tables 11 and 12 
while the combinations Eq. (lll2p are reported in Table 13, for 9 = ^ 2it together 

with the results of 2-parameter fits to the simulated data. Notice the strong scatter of the 
simulated data at the same 0— values. Of course, to compare with the real data, one should 
first take the even combination Eq. (|5Up of the entries in Table 1 which otherwise also contain 
odd-harmonic terms. 

We conclude this brief analysis by emphasizing the importance of Hicks' observation (see 
Sect. 3) concerning the fringe shifts from the session of July 8th. By accepting his interpre- 
tation, the experimental azimuths from the three noon sessions of July 8th, 9th and 11th, 
respectively 9q^^ ~ 357, 285 and 317 degrees, would become Oq"^^ ~ 267, 285 and 317 degrees 
and thus be in rather good agreement with the simulated azimuths reported in Table 13. 
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Table 11: The coefficients 2C(t) Eqs. from a single simulation of 6 rotations in one noon 
session of the Michelson-Morley experiment. The stochastic components of Vx{t) and Vy{t) in 
Eqs. / fgO|) and [6l\) are controlled by the kinematical parameters (V^,a, 7)cmb as explained in 
the text. 



i 


1 


2 


3 


4 


5 


6 


1 


-0.023 


+0.002 


-0.024 


+0.001 


-0.004 


-0.003 


2 


+0.003 


-0.011 


-0.000 


-0.006 


-0.021 


-0.034 


3 


-0.001 


-0.001 


+0.000 


-0.009 


+0.002 


+0.007 


4 


-0.002 


+0.003 


-0.008 


+0.002 


-0.060 


-0.030 


5 


+0.002 


-0.017 


+0.002 


-0.001 


+0.003 


-0.008 


6 


-0.007 


-0.006 


-0.059 


-0.013 


-0.008 


-0.047 


7 


-0.020 


-0.001 


-0.019 


-0.000 


-0.003 


-0.003 


8 


-0.011 


-0.001 


-0.011 


-0.002 


-0.026 


+0.001 


9 


-0.015 


-0.000 


-0.008 


-0.001 


-0.008 


-0.022 


10 


-0.037 


-0.005 


+0.000 


-0.002 


-0.003 


+0.003 


11 


+0.003 


-0.022 


-0.015 


-0.005 


-0.003 


+0.002 


12 


-0.002 


-0.049 


-0.023 


-0.016 


-0.009 


-0.006 


13 


-0.001 


+0.002 


+0.000 


+0.001 


+0.003 


-0.001 


14 


+0.003 


-0.003 


+0.003 


-0.023 


-0.001 


-0.019 


15 


-0.012 


-0.034 


-0.013 


-0.001 


-0.001 


-0.011 


16 


-0.004 


-0.017 


-0.004 


+0.002 


-0.010 


-0.010 
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Table 12: The coefficients 2S(t) Eqs. ^4\) from a single simulation of 6 rotations in one noon 
session of the Michelson-Morley experiment. The stochastic components of Vx{t) and Vy{t) in 
Eqs. \60\) and \61\) are controlled by the kinematical parameters 7)cmb as explained in 

the text. 



i 


1 




2 




3 


4 




5 




6 


1 


+0.011 


+0, 


,001 


+0, 


,011 


+0.000 


+0, 


,003 


+0, 


,001 


2 


+0.004 


-0, 


,003 


-0, 


,001 


-0.008 


+0, 


,007 


-0, 


,009 


3 


+0.003 


-0, 


,005 


+0, 


,000 


-0.007 


+0, 


,007 


+0, 


,000 


4 


+0.001 


-0, 


,007 


-0, 


,003 


-0.001 


+0, 


,001 


+0, 


,016 


5 


+0.002 


-0, 


,007 


+0, 


,000 


+0.001 


-0, 


,000 


-0, 


,003 


6 


-0.004 


+0, 


,010 


+0, 


,001 


-0.001 


+0, 


,006 


+0, 


,009 


7 


-0.017 


+0, 


,002 


+0, 


,005 


-0.000 


-0, 


,003 


+0, 


,002 


8 


+0.011 


-0, 


,002 


-0, 


,012 


-0.001 


+0, 


,010 


-0, 


,000 


9 


+0.011 


+0, 


,001 


+0, 


,008 


+0.011 


+0, 


,005 


-0, 


,011 


10 


+0.016 


+0, 


,005 


-0, 


,001 


+0.001 


-0, 


,002 


+0, 


,001 


11 


+0.000 


-0, 


,013 


+0, 


,015 


+0.005 


-0, 


,002 


-0, 


,001 


12 


-0.001 


-0, 


,022 


+0, 


,004 


-0.003 


-0, 


,000 


+0, 


,002 


13 


-0.001 


-0, 


,001 


+0, 


,001 


+0.000 


-0, 


,004 


+0, 


,001 


14 


-0.001 


+0, 


,002 


+0, 


,002 


+0.018 


+0, 


,001 


+0, 


,009 


15 


-0.012 


+0, 


,018 


-0, 


,002 


+0.001 


-0, 


,004 


+0, 


,007 


16 


+0.002 


+0, 


,008 


+0, 


,001 


+0.002 


-0, 


,005 


-0, 


,006 
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Table 13: The fringe shifts — j-^ Eq. for the single simulation of one noon session of the 
Michelson-Morley experiment reported in Tables 11 and 12. The angular values are defined as 
6 = ^Yg!- 27r. The variance of the averages is about ±0.004 for the amplitude and about ±11° 
for the azimuth. 



i 


1 


2 


3 


4 


5 


6 


average 


1 


-0.023 


+0.002 


-0.024 


+0.001 


-0.004 


-0.003 


-0.009 


2 


+0.005 


-0.010 


-0.001 


-0.010 


-0.010 


-0.031 


-0.010 


3 


+0.003 


-0.005 


+0.000 


-0.007 


+0.007 


+0.000 


-0.000 


4 


+0.002 


-0.007 


+0.004 


-0.002 


+0.043 


+0.033 


+0.012 


5 


-0.002 


+0.017 


-0.002 


+0.001 


-0.003 


+0.008 


+0.003 


6 


+0.008 


-0.003 


+0.041 


+0.010 


+0.002 


+0.027 


+0.014 


7 


+0.017 


-0.002 


-0.005 


+0.000 


+0.003 


-0.002 


+0.002 


8 


-0.016 


+0.001 


+0.001 


-0.001 


-0.025 


+0.001 


-0.007 


9 


-0.015 


-0.000 


-0.008 


-0.001 


-0.008 


-0.022 


-0.009 


10 


-0.015 


-0.000 


-0.000 


-0.001 


-0.003 


+0.003 


-0.003 


11 


+0.000 


-0.013 


+0.015 


+0.005 


-0.002 


-0.001 


+0.001 


12 


+0.000 


+0.019 


+0.020 


+0.009 


+0.006 


+0.006 


+0.010 


13 


+0.001 


-0.002 


-0.000 


-0.001 


-0.003 


+0.001 


-0.001 


14 


-0.001 


+0.001 


-0.004 


+0.004 


-0.000 


+0.007 


+0.001 


15 


+0.012 


-0.018 


+0.002 


-0.001 


+0.004 


-0.007 


-0.001 


16 


-0.005 


-0.018 


-0.003 


-0.000 


-0.004 


-0.003 


-0.005 


Af 
Of 


0.009 
279° 


0.005 
259° 


0.009 
266° 


0.003 
285° 


0.011 

255° 


0.013 

272° 


0.008 
269° 
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